Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



TESTS OF THE RADIAL TREMAINE- WEINBERG METHOD 
Sharon E. Meidt and Richard J. Rand 

Department of Physics and Astronomy, 
University of New Mexico, 800 Yale Blvd Northeast, Albuquerque, NM 87131 

Michael R. Merrifield 

University of Nottingham, School of Physics & Astronomy, University Park, Nottingham, NG7 2RD 

Victor P. Debattista 

Centre For Astrophysics, University of Central Lancashire, Preston, UK PRl 2HE 

AND 

Juntai Shen 

McDonald Observatory, The University of Texas at Austin, 1 University Station, C1402, Austin, TX 78712 

Draft version February 2, 2008 

ABSTRACT 

At the intersection of galactic dynamics, evolution and global structure, issues such as the rela- 
tion between bars and spirals and the persistence of spiral patterns can be addressed through the 
characterization of the angular speeds of the patterns and their possible radial variation. The Radial 
Tremaine- Weinberg (TWR) Method, a generalized version of the Tremaine- Weinberg method for ob- 
servationally determining a single, constant pattern speed, allows the pattern speed to vary arbitrarily 
with radius. Here, we perform tests of the TWR method with regularization on several simulated 
galaxy data sets. The regularization is employed as a means of smoothing intrinsically noisy solutions, 
as well as for testing model solutions of different radial dependence (e.g. constant, linear or quadratic). 
We test these facilities in studies of individual simulations, and demonstrate successful measurement 
of both bar and spiral pattern speeds in a single disk, secondary bar pattern speeds, and spiral winding 
(in the first application of a TW calculation to a spiral simulation). We also explore the major sources 
of error in the calculation and find uncertainty in the major axis position angle most dominant. In 
all cases, the method is able to extract pattern speed solutions where discernible patterns exist to 
within 20% of the known values, suggesting that the TWR method should be a valuable tool in the 
area of galactic dynamics. For utility, we also discuss the caveats in, and compile a prescription for, 
applications to real galaxies. 

Subject headings: galaxies: spiral - galaxies: kinematics and dynamics - galaxies: structure - methods: 
numerical 



1. INTRODUCTION 

One of the prime unresolved issues in the dynam- 
ics and evolution of galaxy disks remains the origin 
and evolution of large-scale bar and spiral structure. 
Though the persistence of grand-design spirals has been 
tied observationally to the presen ce of bars or compan- 
ions (jKormendy fc NormanI [r979( l. virtually nothing is 
known about the actual lifetimes of spiral patterns. Ad- 
ditionally, despite indications that t he relation between 
bar an d spiral pattern speeds (which ISellwood fc Sparki 
(jl988) first argued may not be equal) may be im- 
portant for understanding the role of bars in angular 
momentum transfer during secu lar disk evolution (e.g. 
Debattista fc SellwoodI [19981 and iDebattista fc SellwoodI 
2000( 1 ■ there are as yet unanswered questions about the 
connections between multiple patterns in different radial 
zones. While mode-coupling between patterns, which 
allows efficient outward angular momentum tran sfer in 
disks (ISvgnet et al.lll999l : iMasset fc Taggeilll997f ) seems 
a most promising link, in 2D N-body simulat i ons w ith a 
dissipative gas component iRautiainen fc Said (fl999) find 
evidence for spiral structure in the absence of a bar, bar- 
spiral mode coupling, spiral-spiral mode coupling, and 



multiple pattern speeds without mode coupling. 

Clearly, to address questions about the persistence of 
spiral patterns and the relation between bars and spirals 
requires not only determination of the pattern speed but 
how it varies with radius; only with accurate measure- 
ment of bar and spiral or inner and outer spirals pattern 
speeds in the same galaxy can we confirm whether spiral 
structure is steady or winding, whether bars and spiral 
pattern speeds are equal or are unrelated, whether mode- 
coupling exists, and the domain and number of patterns 
that can be sustained in a disk. 

Because they are not directly accessible through ob- 
servation, pattern speeds are often determined with in- 
direct means such as the identifi cation of predicted be- 
havior at resonance r adii (e.g. lElmegreen et al.l Il989t 
lElmegreen et al.l 119961 ) or kinematic and morphological 
comparisons of simulated and observed structur e (e.g. 
IRautiainen et al.ll2005l: iGarcia-Burillo et al.l ll993'). It is 
also clearly desirable to employ methods for estimating 
pattern speeds which do not rely on theoretical mod- 
els or simulation. Many other pattern speed determina- 
tions have therefore centered on the use of the model- 
independent method of Tremaine & Weinberg (1984; 
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hereafter TW) which presents a rigorous derivation for 
the pattern speed Hp based on the requirement of conti- 
nuity using observationally accessible quantities. The de- 
termination of flp involves surface density- weighted posi- 
tion and velocity line integrals parallel to the galaxy ma- 
jor axis under several essential assumptions. Specifically, 
the method requires that the disk of the galaxy is flat 
(unwarped) and contains a single, well-defined rigidly- 
rotating pattern; that the surface density of a kinematic 
tracer of a disk component, which must obey continuity, 
becomes negligibly small at some radius and all azimuths 
within the map boundary (thereby critically yielding con- 
verged integrals); and that the relation between the emis- 
sion from this component and its surface density is lin- 
ear, or if not, suspected deviations from linearity can be 
modeled. 

In that information from all sampled radii is associated 
with a single, constant pattern speed, the TW calcu- 
lation poses a challenge for extracting multiple distinct 
or radially varying pattern speeds. Non-axisymmetric 
structures beyond a dominant pattern such as a bar in 
the disk of a galaxy will interfere with the measure- 
ment of rip of the bar; when a non-axisymmetric disk 
can be decomposed into two components with different 
pattern speeds, then the TW estimate is a luminosity- 
and asymmetry-weighted average of t he two patterns 
(jPebattista. Gerhard fc Sevensteii[200l . The TW esti- 
mate for the secondary pattern speed in the inner disk 
of NGC 6946 observed in Hq., for example, is estimated 
to be limited to an uncertainty of as much as 50% given 
the primary patte rn's contribution to the TW integrals 
(jFathi et al.ll2007f) . 

Such issues notwithstanding, there have been several 
recent adaptations of the TW method to positive ef- 
fect. Applications of the TW method to SBO galaxies 
using stellar light as a tracer extend the limits of inte- 
gration in the TW calculation to just past the end of 
the bar in order to minimize contributions to the TW in- 
tegrals from non-axisymmetric (and several magnitudes 
dimmer) features beyond the bar; in such cases, inte- 
grating past the structure of interest is found to be suffi- 
cient for achieving con verged integrals. The bar pa ttern 
speeds in N GC 7079 (iDebattista fc Wilhamsll200l and 
NGC 1023 (jPebattista et all I2002| ). for example, have 
both been successfully measured in this way. 

When there exist more than one pattern in distinct 
radial zones, however, arguments about the convergence 
of the TW integrals are less straight-forward. To mea- 
sure the secondary bar pattern spe e d in NGC 2950, 
iCorsini et all (|2003f) and lMacieiewskil (|2006f) explore de- 
coupling the inner secondary and outer primary bar pat- 
tern speeds by associating each component with unique 
surface brightness contributions in the TW calculation. 
The TW integrals are modified by the presence of the 
inner pattern based on an assumption about how and 
where the two patterns decouple (to first order). This 
analysis has confirmed the existence of, if not measured, 
a unique secondary bar pattern speed possibly indicating 
counter-rotation with respect to the primary bar. 

At their best-improving the accuracy of pattern spee d 
estimates to about 20% (jCerssen fc Debattistal l2007( )- 
these kinds of adaptations of the TW method for mea- 
surement of single or multiple patterns still require as- 
sumptions about bar extent based on morphological or 



kinematic signatures. To separate the observed surface 
brightness in NGC 2950 in t o seco ndary and primary bars 
components, iMacieiewskil ()2006[ ) must assume that the 
patterns indeed decouple at the inner-bar end, or that 
the outer pattern is axially symmetric at least within 
the inner's extent. There, the transition between the 
two is inferred from the location of a plateau in the TW 
surface brightness-weighted position integral as the lim- 
its of integration are extended from zero. Perhaps more 
critically, as investigated in § 14.21 the direct association 
of bar length measured in this manner (or perhaps oth- 
ers) with pattern extent may introduce error into TW 
calculations. 

Furthermore, adaptations of the TW method based on 
this type of identification are likely to be inapplicable for 
spiral pattern estimation. In spiral galaxies, not only can 
identifying transitions between patterns be less clear, but 
non-axisymmetric motio ns are significantly smaller (e.g. 
[Roberts fc Stewartlll987f ) than in bars (at least those typ- 
ically analyzed with the TW method). Indeed, the TW 
method has yet to be tested on a simulated spiral. 

A recent modification to the TW calculation (Merri- 
field. Rand fc Meidt 2006; hereafter MRM) in which ftp 
is allowed radial variation promises to be an invaluable 
resource for tests of long-lived density wave theories and 
for understanding the connection, if any, between bar 
and spiral pattern speeds. Like the TW method, the 
so-called radial TW (TWR) method uses measurements 
of observables to extract radially varying pattern speeds. 
As first applied using the BIMA Survey of Nearby Galax- 
ies CO observations of the grand-design Sb galaxy NGC 
1068 (MRM), the TWR method returned a spiral pattern 
speed solution that declines with radius, allowing a wind- 
ing time for the pattern to be estimated (e.g. MRM). 

As described in MRM, the nature of the discretized 
calculation presents numerical solutions that are highly 
susceptible to fluctuations as a result of compounded 
noise in the data. In this paper 12. 2p . we develop the 
TWR method with regularization as a means of smooth- 
ing intrinsically noisy solutions as well as testing model 
solutions of different radial dependence (described in § 
12. 3p . These and other commendations notwithstanding, 
with regularization one may risk introducing an unreal- 
istic prejudice to TWR solutions. In § 13.11 we address a 
means of identifying when this is likely to occur, and in § 
I3.2l describe a scheme for minimizing such regularization- 
induced bias. 

Using evidence that arises from these considerations, 
as well as other a priori information, theoretically and 
observationally motivated models for f2p(r) can be devel- 
oped which constrain the number and extent of patterns 
present in a disk. Once solutions for these models are cal- 
culated, the goodness of each must be assessed in order 
to identify the best-fit solution. In § 13.31 outline the 
criteria with which the models are judged and describe 
our concept of error evaluation in the final solutions. 

Beginning in § 14.11 we analyze three simulated galax- 
ies with known pattern speeds (a barred spiral, a slowly 
winding spiral, and a double barred spiral, marking the 
first application of the TW method to simulated spiral 
patterns) in order to develop a general stratagem for ap- 
plication to real galaxies. As applied to these simula- 
tions, we find that the TWR method is able to extract 
multiple pattern speeds with accuracies on the order of 
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(and, as we will see in some cases, better than) the tra- 
ditional TW method. With the barred spiral simulation 
in § 14.21 we show that TWR bar pattern speed measure- 
ment presents an improvement over traditional TW bar 
estimates, particularly when there is evidence of a sig- 
nificant contribution from the spiral pattern to the TW 
integrals. We find that the regularized TWR method 
can recover information from both patterns effectively 
by identifying and treating the bar-to-spiral transition 
radius (which the TW values themselves may not indi- 
cate) as a free parameter in the calculation. In § I4.2.1i 
we analyze the results of the method in detail, particu- 
larly with regard to morphological limitations. We com- 
pare our TWR results with TW estimates in § 14.2.21 and 
examine the infiuence of systematic errors due to the as- 
sumed disk position angle a nd inclination (sh own to be 
crucial for TW estimates in iDebattistal (|2003[ )) on both 
in § 14.2.31 We also explore the reliability of a fixed pa- 
rameterization for the bar-to-spiral transition radius. (§ 

11231). 

In the last third of the paper, we investigate the 
prospects for extracting spiral pattern speed solutions 
that are winding in nature (§ 14. 3p . marking the first 
application of the TW method to a spiral simulation, 
and in § 14.41 we address the use of the TWR method for 
the purposes of parameterizing an independently rotat- 
ing nuclear bar in a double barred simulation. There, 
the techniques we employ for decoupling and extracting 
measurements of the pattern speeds of both the primary 
and secondary bar components may present an interest- 
ing corollary to recent attempts with the TW method 
to measure secondary bar pattern speeds in the presence 
of a strong primary bar pattern. We note here that the 
TWR method is a gene ralized version of the procedure 
used on NGC 2950 (see lCorsini et al.ll2003t iMacieiewskil 
[2006); separating the surface brightness into two com- 
ponents can be thought of as the coarsest version of the 
discretization that is the back-bone of numerical TWR 
solutions. 

Based on our experience with these simulations, we 
conclude with comments on the applicability of the 
method to observations of real galaxies in § 15.11 where 
we also outline a general prescription for using the TWR 
calculation with regularization. 

2. THE TWR METHOD WITH REGULARIZATION 
2.1. The Radial Tremaine- Weinberg Method 

By proceeding under the aforementioned assumptions 
of Tremaine & Weinberg (1984), but allowing that flp 
may possess spatial variation in the radial direction 
whereby the surface density of the chosen tracer can be 
written E(a:,y,i) — ^{r,<p — r2p(r)t), with appropriate 
mathematical generalizations, the derivation and mea- 
surement of the pattern speed flp{r) can be made from 
observable intensities and kinematics of a chosen tracer. 

Following the derivation given in MRM, integrating the 
continuity equation obeyed by the tracer (with the re- 
placement dYi/dt = ~Q,p{r)dYi/ d4)) 



dx 



dY^Vy 
dy 







(1) 



over X and y (thereby eliminating the unobservable Vx 
and the spatial derivative d/dy), and changing from 



Cartesian to polar coordinates yields 

l^oc p-K—arcsin{y /r) poo 

/ V,p{r)-^rdrd(f> + / 'Evydx = 0. 

I r—y J cj)—arcsinijj /r) J —oo 

(2) 

A final integration with respect to results in a Volterra 
integral equation of the first kind for ^p{r): 

/oo 
Y^Vydx 
-OO 

(3) 

where x' [r^ y) ~ — y^ ■ 

Note that with constant ^p in equation ([3]), we arrive 
at the regular TW result 



f2p / Y^xdx 



T,Vydx, 



(4) 



which, with normali zation by / T}dx (e.g. 
iMerrifield k Kuiiken|[l995[) . leads to 



< V > 

< X > 



(5) 



where <v>— J Yvydx/ J Yidx and <a;>= / T^xdx/ J Y^dx. 

For a galaxy projected onto the sky plane with inclina- 
tion a (so to distinguish from index i), both the kernel 
on the left and the integral on the right of equation ^ 
are observationally determined quantities with x = Xots , 
y = Hobs/ cos a, and Vy = Vobs/ sin a, where Xobs and yobs 
are the coordinates in the plane of the sky along the ma- 
jor and minor axes, respectively, and Vobs is the observed 
line-of-sight velocity. Solutions can be extracted numeri- 
cally by replacing the integral on the left with a discrete 
quadrature for different values of y = yi and r = rj (see 
Figure [J) whereby equation ([3]) is converted to 



^r,>y,K{yi,rj)Qp{rj) = b{yi) 
or to a matrix equation of the form: 



(6) 



(7) 



with K an upper triangular N x N square matrix. Note 
that numerical quadratures on either side of the galaxy (y 
<0 or or y >0) occur independently, providing two mea- 
sures of flp{r). Furthermore, as governed by the infor- 
mation available, the slices which delimit the quadrature 
on a single side need not be uniformly spaced. In this 
case, solutions inherit a variable bin width Ar. Also, the 
calculation allows for no azimuthal dependence for the 
pattern speed, which we assume throughout. 

The size of K depends on the desired coarseness or 
fineness of the quadrature; the separation between slices 
at positions yi (limited by either the resolution or the 
sampling of the data) translates into a radial bin width 
(modulo cos a) via equation ([7|). The quadrature, per- 
haps more critically, depends on the limits of integra- 
tion in equation ([3]). These limits iX^ax should be 
chosen based on where the integrals have converged. 
While in the case of a single bar pattern integrating past 
the structure of interest is often suitable, as shown in 
iZimmer. Rand fc McGrawl (|2004) (Figures 9-11), in the 
presence of strong, extended asymmetry TW values are 
highly dependent on the extent of integration along each 
slice i. In cases where multiple patterns exist in a single 
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X (kpc) 



Fig. 1. — Illustration of a y>0 quadrature for a galaxy viewed 
from the negative z side with a tilt around the y axis by a=45°. 
The horizontal lines, or slices, at positions yi are spaced at 
Aj/=1.54, or Aj/o6s=l-09, and represent integration between the 

limits ± ^J /Jq + y^ where _Ro=10.8 is the maximum radial extent 

of the quadrature. Each slice is carved into elements of width Ar 
whereby all the elements with the same shade of grey represent a 
single radial bin rj. 

disk, then, it is equally favorable (and hopefully suffi- 
cient) to extend all integrals to the edge of the surface 
brightness distribution. 

Meeting the requirement of integral convergence in this 
manner as applied to the TWR calculation determines 
the location of the last radial bin jmax associated with 
elements Kij^^^ along each slice. For a given radial bin 
width, with the requirement that jmax equals N we are 
presented with the size of K as well as the outermost slice 
position, since jmax must also equal imax- One should 
check to see that Knn, the last entry in K associated 
with the outermost slice, is associated with a fully con- 
verged / Yjxdx (which achieves convergence at least by 
the map boundary). 

Since K is an upper triangular matrix, the fJj can be 
solved for via simple back-substitution. In this way, solu- 
tions are generated from the outermost to the innermost 
radius (from light gray to black in Figure [T]) according to 



(8) 



with k > i (from eq. [7]). Such solutions, however, are 
especially susceptible to wildly oscillatory behavior as 
errors from large radii propagate inward. These com- 
pounded errors arise from uncertainties in the velocities 
and intensities (which get translated into the Kij) and 
can be particularly severe since the outermost bins often 
cover the lowest signal-to-noise regions in the data. The 
errors thus introduced can be systematic. 

Figure [2] shows an example of (unregularized) TWR 
solutions for two different binnings using data from the 
barred spiral simulation of § 14.21 By increasing the bin 
width the largest oscillations are reduced, but (as will be- 
come clearer) even tripling the bin width will not neces- 
sarily limit the propagation of noise to the level required 
for the extraction of realistic solutions. (A more thor- 
ough discussion of the implications for TWR solutions 



will be deferred until the beginning of § 14.21 ) 

Though the initial application of the TWR method on 
a real galaxy, namely NGC 1068 (MRM), showed little 
of the oscillatory behavior common to noisy, discretized 
Volterra-type solutions (outside kpc, anyway), the 

solutions were generated over relatively few bins (only 
five at most, over the region r—1.5 to 2.8 kpc). In gen- 
eral, while large bin widths can often minimize the prop- 
agation of noise in the calculation, they can be expected 
to compromise solutions, as will be discussed in the sec- 
tions to follow; naturally, the smaller the quadrature el- 
ement, the more accurately the true radial variation of 
the pattern speed can be ascertained. This is particularly 
critical as applied to disks sustaining multiple patterns 
where, as described later, a certain degree of radial pre- 
cision is required for accurate separation of the pattern 
speeds. 

Barring a large, limiting resolution, one should expect 
to be able to perform a sufficiently smooth quadrature 
wherein the number of elements in K becomes large. 
Since more elements in K (and more bins over which to 
generate solutions) result in intrinsically noisy behavior- 
an effect most pronounced in the inner-most bins-gaining 
a finer, more accurate quadrature often means forfeiting 
control of the solution. By combining regularization with 
the TWR calculation to force a smooth solution, how- 
ever, one can counter this effect while maintaining the 
required precision. 

Regularization also serves to alleviate the impact of 
non-global features that are most likely not included in 
the overall pattern (and which can singularly introduce 
large errors into the integrals). Since rapid fluctuations 
in f2p(r) are penalized, discrepant points need not be 
avoided or ignored (as demanded in performing the TWR 
calculation on NGC 1068 in MRM). While one may also 
use the alternative, which would be to fit models of Qp 
directly to equation ([7]) and perform a grid-search to find 
the best model form and coefficients, we pursue regular- 
ization here, its speed making it preferred. 

2.2. Regularization 

Our procedure entails the following. As a modifica- 
tion to the estimator minimized by solutions flj of 
equation ([7]), namely 



iKijftj - bi 



(9) 



with implicit sum over i (and j) and errors ai repre- 
senting the measurement error of the i*'' data point bi, 
we introduce a regularizing operator, or smoothing func- 
tional S, containing a priori inforr nation in the ma nner 
of Tikhonov-Mil ler regularization (jTikhonov fc ArseninI 
Il997t lMilleii[T97Q) in which (in matrix form) solutions Jl 
minimize 

|K-rj-bp-|-Af2-S-r2. (10) 

Here, the elements of K and b are Kij/ai and bi/ui, 
respectively, and the role of A-controUing the relative 
amount of minimization on the left to entropy maxi- 
mization on the right-is explicit. 

Reduced to a linear set of normal equations, this mini- 
mization returns smoothed solutions according to a mod- 
ified version of equation ([7]): 

(K^ •!<: + AS) • n = • b (11) 
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Fig. 2. — Pl ots o f (unregularized) TWR solutions for two different binnings of data from the barred spiral simulation in 
SA=-45° (see 14. II for orientation convention). The left (right) panel shows the solution generated using Ar=0.3 (0.9) bins. 
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Note that the regularizing functional, not necessarily up- 
per triangular, introduces an anticipatory quality to solu- 
tions flj whereby all bins at the same radius are coupled. 
Furthermore, solving for components Clj no longer only 
involves a procedure like back-substitution, but requires 
rather an L-U decomposition (for instance) as well. 

2.3. The Smoothing Operator 

The real power in applying regularization to TWR cal- 
culations is in the freedom to choose how the smooth- 
ness of solutions is achieved. For the purposes of distin- 
guishing between different possible radial dependences 
for r2(r), we choose S to reflect a priori assumptions 
based on simple expectations from theory and observa- 
tion. Model solutions, then, each incorporating its own 
S, represent smoothed, testable realizations of the pat- 
tern speed. These we restrict to simple forms in order 
to minimize the additional amount of information to be 
extracted from the data relative to the traditional TW 
method. 

For polynomial solutions, we consider only constant, 
linear and quadratic radial dependence. The elements of 
the smoothing S are associated with the minimization of 
the rt*'* derivative of il{r) for each polynomial solution 
of order n. For instance, for linear solutions this entails 
minimizing 



ft s n = 
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One may also choose a form for S that identifies two 
or more distinct regions of independent radial behavior 
by invoking a step-function model. For the case of a 
barred spiral with a constant bar and quadratic spiral. 



for instance, this corresponds to minimizing 



t-i 

E 

7V-3 

El 

n—t 



n+l\ 



for n<t 



Q„ + Sfln+l - Sfln+2 + ^n+sl'^ for n><. 

(14) 



The elements of S with n < t reflect the a priori assump- 
tion that the bar pattern speed is constant, while those 
for n > t associate a quadratically varying pattern speed 
with the spiral. The index t, a free parameter, locates 
the radial bin where the transition between the two pat- 
terns occurs. Obviously, the number of available bins 
constrains the order of the polynomial in a given radial 
zone. 

Once we have chosen S, we initially choose A to reflect 
comparable amounts of minimization and regulariza- 
tion by letting A = Aq = Tr |k'^ • k| /Tr{S}. Since 

we are in the business of generating solutions based on 
particular models, A is modified to arrive at the regu- 
larization required to return solutions of a given type. 
This modification generally consists of an increase in A 
over Ao . Consider how the regularizing parameter A reg- 
ulates the degree of smoothness of the solution to the 
weight placed on the data: with A=0, equation (fTO|) cor- 
responds to minimization (and becomes an unbiased 
estimator with the smallest variance), however yielding 
highly oscillatory solutions, while A— >cx) corresponds to 
a maximally smooth estimator with non- vanishing vari- 
ance. 

Fitting data sets with different spatial coverage will 
change the effect of A on the solution (e.g. larger bins 
require less regularization) . The most appropriate choice 
for A (and S) should be made on a galaxy- by-galaxy 
basis, according to the quality of information to be ex- 
tracted from observations. 



3. 
3.1. 



OTHER CONSIDERATIONS 



Regularization-induced Bias 

By imposing assumptions about the smoothness of the 
pattern speed, regularization inevitably introduces com- 
plications for extracting realistic solutions. To under- 
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stand how these arise, consider the sohition for a barred 
spiral galaxy. The nature of the calculation (from out to 
in) has implications for the accuracy of the bar estimates, 
in particular. Not only do the ilj for bins covering the 
bar rely on the greatest number of matrix elements (and 
errors therein), but the bar estimate depends critically on 
the solution for the spiral and all outer bins via equation 
([8]). Consequently, merely requiring the pattern speed in 
the outer bins to be constant with regularization out to 
the edge of the surface density (for instance)-effectively 
removing fluctuations that might better fit the data-will 
have consequences for the bar solutions. So while the 
regularization is particularly fast and effective for tests 
for the radial behavior of patterns, it can also hinder the 
realization of accurate solutions. 

For the simulations studied here, the risk of 
regularization-induced bias is inherited from the adopted 
quadrature. Recall our requirement that all slices cover 
the full extent of the 'emission'-so as to insure all inte- 
grals be fully converged-and relatedly, that the last ma- 
trix element governs the outermost slice position. In the 
barred spiral simulation of § I4.2i for example, such an 
extensive quadrature presents us with outermost slices 
that pass through a region where there is simply no dis- 
cernible pattern (as indicated by the surface brightness 
distribution and its Fourier decomposition; see next sec- 
tion) . While these slices themselves do not provide direct 
estimates of the patterns of interest, the corresponding 
bin values are necessary for calculating the bar and spiral 
solutions. Moreover, the quality of these solutions will 
be intimately related to the treatment of the outermost 
bins. We therefore find that identifying, and reducing 
the influence of, the compromised zone by not enforcing 
regularization on these bins to be essential for accurate 
pattern speed measurement. 

3.2. Fourier Diagnostics 

With the above concerns in mind, we have tested and 
used the following scheme. Given slices that pass through 
an outer region which either contains little information 
from a strong pattern, is suspected of sustaining multi- 
ple patterns, or displays only faint emission, we choose 
in such a case to let the values in the outermost bins 
be calculated without regularization with the restriction, 
only, that they minimize the x^- Once a particular bin- 
at Tc, the cut radius-has been reached, regularization is 
imposed with all remaining inward bins generated ac- 
cordingly. 

In our procedure this corresponds to using an S in- 
dexed by the cut bin c that is lowest-block zero. And 
the 'cut' radius identifies the location in the disk where 
the outermost discernible pattern ends. Note that this 
is in contrast to altogether ignoring the outer portion of 
the disk. We prefer this procedure for two reasons: 1) as 
qualified later in § 14.2.21 when each integral is truncated 
within the disk, the quadrature is at greater risk of ig- 
noring information critical for characterizing the patterns 
uniformly throughout the disk and 2) like the transition 
radius r^, we can easily incorporate rc as a free, though 
restricted, parameter in our models. 

In practice, extracting the bar and spiral pattern 
speeds for the simulated barred spiral in § 14.21 involves 
generating a group of solutions with various bar-to- 
spiral pattern transition radii for a given 'cut' location. 



Throughout the analysis, we choose the cut radius to 
reflect a priori knowledge of the outermost measurable 
pattern's termination radius estimated from the surface 
brightness and its Fourier decomposition. When referred 
to, the power in each Fourier component, or mode m, is 
given by the norm of the complex Fourier amplitude 

N 

I, = e™^" (15) 

n=l 

where 0„ is the angular coordinate of each of the N par- 
ticles at each measured radius. 

For the barred spiral in § 14.21 for example, we com- 
bine evidence from the surface density-where beyond the 
bar there is enhanced spiral surface density from only a 
limited radial zone-with the Fourier spectrum to iden- 
tify a region in the disk outside the spiral which is sus- 
ceptible to regularization-induced bias. Specifically, the 
Fourier power spectrum (Figure |4|) shows that at r~5.0 
the second clean hump in m=2 power decreases to almost 
zero, marking the end of the spiral. Past this radius, 
the (strong) m—2 component (between r~6. 0-8.0) is not 
associated with visible spiral structure (see Figure [3]). 
Reckoning this outer zone to be incompatible with a sim- 
ple pattern speed model, then, we consider only the in- 
ward bar and primary spiral pattern spe eds to be measur- 
able w ith regularization. Figure 19 of ' Debattista et all 
(|2006| ) for this simulation (Figure [5] in this paper) con- 
firms this; not only does the spiral pattern terminate 
at rc~ 5.0, but beyond this radius the pattern speed is 
multi-valued (this, of course, would be indiscernible in a 
real galaxy). In this case, imposing form with regular- 
ization on bins of suspect quality and behavior outside 
the spiral will likely impair the solution of interest. We 
therefore restrict the cut bin for the barred spiral simula- 
tion to 4.5<rc<6.0, representative of where the primary 
spiral pattern terminates in the disk. As mentioned in 
§ 14.21 this step is substantiated by our finding that a 
cut radius of rc=4.8 is one of several minima given a 
range of possible cut radii. And furthermore, solutions 
generated in this manner are judged to overall provide a 
considerably better fit to the data than solutions where 
regularization is imposed out to the edge of the surface 
brightness (according to the scheme described in the next 
section) . 

3.3. Weighting Schemes and Goodness-of-Fit 

Given the data, we simultaneously generate model so- 
lutions with different radial dependences for direct com- 
parison for each side (y>0 or y<0) independently. We 
then average the two like-model solutions together to 
construct a global solution. (In those instances when 
model solutions include a cut bin, the averaging occurs 
over the regularized part of the solutions only, in order 
to maintain the 'unbiased' quality of the unregularized 
part of the solutions for each side of the galaxy.) Note 
that while the assumption that the patterns are indeed 
global is not overly inspired for the simulated galaxies 
studied here, putting this into practice on a real galaxy 
requires that the assumed galaxy kinematic parameters 
are accurate and that such symmetry exists. 

With each global model solution we generate a com- 
plete set of <v> using equation ([7]). We then judge 
each model through the xl ix^ per degree of freedom) 
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goodness-of-fit estimator of the reproduced to actual 
<v>i= bi/{J 'Edx)i given measurement errors crf^^, in 
keeping with the standard TW analysis. Note that with 
this choice the ai in equation ([9]) are related to the er- 
rors af^-^ by (/ T,dx)i. That is, the calculation fits to 
the bi given errors Oi while our estimator considers 
the differences from <v>i given errors <jf^^ . 

While for a real galaxy inaccuracies in the assumed 
position angle (PA) have the largest potential for intro- 
ducing errors into <t)>i, we prefer that the measurement 
errors a^^^ reflect random noise in the data, only. (Sys- 
tematic errors prove more practically assessed through 
direct tests of the sensitivity of the results to departures 
from the nominal values for the PA or inclination, for 
instance.) For the simulations studied here, then, we ob- 
tain errors (7^^^ under the assumption that the inverse 
mirror image of each <v>i on one side of the galaxy 
should be the same as on the other side (i.e. the patterns 
are symmetric). We then assign a global error cr<"> to 

each slice where cr<^> = (x;«fi '^t"^) /2N (and N is the 
number of bins/slices used in the TWR calculation on a 
single side). 

In practice, the simplest weighting schemes are ei- 
ther uniform weighting for all slices or weighting by the 
intensity, which should give more weight to slices where 
the signal is strongest. We have chosen the former since 
we are interested in fip over a broad radial range, and 
prefer that our result is not dominated by just the slices 
with the highest signals (which can vary dramatically). 
Furthermore, the choice of assigning an identical error 
cr^"^ to each slice carries with it an implicit weighting 
scheme for equation pO)) . For an exponential surface 
brightness profile, for example, slices on either side of 
the galaxy will have progressively smaller J Sdx as \yi\ 
increases. This corresponds to errors cr,; proportional to 
J Yidx^ then, that grow larger from out to in. Since in 
most cases the uniform weighting scheme will be in ac- 
tuality most restrictive of the outer bins-the goodness of 
which will affect the solution inward-this choice is par- 
ticularly well suited for the TWR calculation. 

Given our choice of weighting scheme, there are two 
important considerations which demand that we calcu- 
late the over all slices. First, the inner-most bins 
contribute to the <v> in only the inner most slices and 
hence contribute relatively negligibly to the x^, despite 
possibly larger weights Kij (reflective in part of a sur- 
face brightness that is centrally peaked, say) than bins 
at larger radius. This is especially true when an inner 
pattern appears only over a small fraction of the total 
bins, and the reproduced <v> of even those slices that 
pass directly through the inner pattern still rely (per- 
haps predominantly) on the solution out to the largest 
radial bin. However, since the goodness of the inner bins 
is directly related to the goodness of the bins at larger 
radii, by considering all slices we more effectively judge 
the whole solution. 

Secondly, in cases when some number of outermost 
bins are calculated without regularization, it is critical 
to account for the (largely positive) effect that these bins 
have on the solution inward, especially from model to 
model. In our current scheme, the values of the 'un- 
regularized' bins are not quite identical to those in the 
completely unregularized solution since the quantity that 



they minimize still includes participation from non-zero 
elements in the smoothing functional S and what is cur- 
rently a model-dependent A. That is, slight variation 
in the values of the 'unregularized' bins from model to 
model is apparent, and we cannot ignore the minor dif- 
ferences this introduces to the regularized part of the so- 
lution. Though this is a minor effect, by considering the 
<v> of all slices in the we prevent solution-preference 
based on the 'unregularized' bins (which act essentially 
as stand-ins) from being introduced. 

The first consideration above also prevents us from cal- 
culating judicious error bars on the solutions according 
to the variation of individual model parameters over a 
typical confidence interval; in practice, the value of 
an inner pattern with minimal radial extent can change 
considerably with little effect on the x^. Indeed, we find 
that the errors generated according to such a prescrip- 
tion are unrepresentative of the goodness of the solutions 
as returned by the calculation. 

In § 14.2.11 we describe the dependence of both inner 
and outer speeds on the location assigned to the tran- 
sition between the two. An obvious progression for fu- 
ture applications of the TWR method would be an explo- 
ration of the covariance of what we consider here, to first 
order, 'free parameters', especially for the purposes of 
improved error estimation. Presently, however, we con- 
struct error bars for model solutions by considering the 
range of parameters in the best solutions at different as- 
sumed projections. As will pertain to the sections that 
follow, by considering an overall solution in this manner 
we can fairly account for the uncertainty introduced for 
real galaxies by the reality that each can be sampled at 
only a single PA. 

4. TESTS OF THE TWR METHOD ON SIMULATIONS 

In order to establish guidelines for applying the reg- 
ularized calculation to observations of real galaxies, we 
next perform tests of the method on simulations with 
known pattern speeds. Each case invokes unique models 
for r2p(r) which we motivate and discuss in detail. The 
procedure for engaging the method with maximum accu- 
racy then follows from careful examination of the qual- 
ity of solutions given the available information. Though 
quite detailed, these individual studies together constrain 
general scenarios and practices to extrapolate onto ob- 
servations of similarly structured, real galaxies. 

4.1. N-body Systems 

We use three simulations in this study. The first, which 
we refer to here as simulation I, constitutes a barred 
galaxy with spiral stru cture. Originally presented in 
iDebattista et al.l ()2006f ) (hereafter D06) where it IS re- 
ferred to as run L2.tl2, it consists of a live disk immersed 
in a rigid halo. A complete description of the model pa- 
rameters are contained in D06. 

The second (spiral) simulation, which we refer to here 
as simulation II, is unpublished. It was designed with 
the main aim of generating stro ng spiral structur e using 
the groove mode mechanisn i of ISellwood fc LinI (|1989l ) 
and 'S ellwood fc Kahnl (|1991| ) in which dynamical insta- 
bility develops from a 'groove', or narrow feature, in the 
phase-space density at a particular angular momentum. 
A trailing spiral wave is generated, and at the Lindblad 
Resonances of the wave, further grooves develop such 



8 



o 




o 



-10 10 

X 

Fig. 3. — Face-on display of the barred spiral simulation's surface 
brightness distribution projected with a -30° rotation about the z 
axis. For reference, the alignment of the TWR quadrature for a 
frame at this orientation is designated SA=-|-60°. 

that the instability is recurrent. Like simulation I, it 
consists of a rigid halo and live disk, but also includes a 
live bulge component. The bulge constitutes 25% of the 
baryonic mass and is sufficiently concentrated that a bar 
is very slow in forming. The disk has Toomre-Q — 1.2; in 
order that a strong spiral was seeded, ~ 6% of disk parti- 
cles in a narrow angular momentum range were removed 
leaving 4x10^-169480 particles including the bulge. The 
result, as can be seen in Figure [TTl is the formation of a 
strong but transient spiral. 

The last simulation in this paper, simulation III, is 
a double-barred ga l axy g enerated using the method of 
iDebattista fc ShenI (j2007[ ). This high- resolution simula- 
tion consists of live disk and bulge components in a rigid 
halo potential. The model has ~ 4.8 million equal mass 
particles, with ~ 4 million in the disk and ^ 0.8 million 
in the bulge such that the bulge has mass Mf, — 0.2Md, 
where Md is the disk mass. The initial Toomre-Q of 
the disk is ~ 2. The formation of the secondary bar is 
induced by making the bulge rotate (to mimic a pseu- 
dobulge). More details of the simulation can be found in 
IShen fc Debattistal ()2007[ ) where it is referred to as run 
D. 

As in D06, all lengths and velocities are here presented 
in natural units. We analyze a snapshot of each simula- 
tion at a single time step with the disk in the xy plane. 
By rotating the system about the z axis we assign a line- 
of-sight direction to establish the kinematical major axis. 
Another rotation about the x axis gives the system an 
inclination a (chosen throughout at q;=45°, unless other- 
wise specified). The snapshot is then projected onto the 
sky-plane where Xobs—x and yobs=y cos a. For a given 
slice spacing, the slices along which the calculations oc- 
cur are aligned perpendicular to the line-of-sight direc- 
tion (parallel to the kinematical major axis). The orien- 
tation of these slices, which is identical to the disk PA in 
a real observation, is designated uniquely the slice angle 
(SA) in the studies that follow. 

4.2. Simulation I: Barred Spiral Galaxy 
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Fig. 4. — Fourier power spectrum of the barred spiral simulation's 
surface brightness distribution shown in Figure |3] Modes up to 
m = 4 are plotted as a function of radius with lines for m = 1 in 
dot, m = 2 in solid, m = 3 in dash-dot-dot, and m = 4 in dash. 
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Fig. 5. — Contours of the barred spiral simulation's m=2 Fourier 
mode showing a bar pattern speed r2i,=0.29, a bar-to-spiral transi- 
tion of rt~2.5, and a dominant spiral pattern speed f2s=0.18 out 
to r'^4.5-5.5 beyond which multiple spiral modes exist. 

The bar and spiral structure in this simulation, first 
presented in D06, is featured out to r~5.0, clear in the 
surface density (Figure [3]) and its Fourier decomposition 
(Figure [4]). Beyond r~5.0, the Fourier decomposition 
indicates the possible presence of a third pattern. With 
step- models for flp{r), then, we might reasonably extract 
pattern speeds for three distinct structures. However, the 
m=2 mode between 5.0<r<8.0 is not associated with a 
strong surface density enhancement. And as remarked 
upon in i) l3.21 Figure 19 in DOG reproduced here in Figure 
[5] shows that the pattern in this radial zone is maintained 
by multiple distinct pattern speeds. (Note that a real 
galaxy would not be disposed to the analysis provided 
with this type of plot. It is available for this simulation, 
only, and we include it here for the sake of compari- 
son). In a clear account of regularization- induced bias, 
test solutions based on a three-pattern speed model have 
considerably larger than those parameterizing a bar 
and single spiral. We therefore reject models with a third 
pattern speed and use our 'cut' scheme where solutions 
are generated without regularization up to a cut radius 

which parameterizes the end of the primary spiral pat- 
tern. 

The actual pattern speeds of the bar and spiral struc- 
ture to be reproduced by our solutions are fli,— 0.29 for 
the bar, and a constant f7s=0.18 for the spiral, as es- 
timated from Figure [5] (from which we also estimate a 
bar-to-spiral pattern transition radius r~2.5). In our 
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Fig. 6. — The best- fit solution and error bars from the TWR 
method apphed to the barred spiral simulation averaged over six 
SAs. The bar C5=0.31±0.02 and spiral r2s=0.17±0.03 are shown 
as solid lines with dashed errors. Errors in the bar- to-spiral transi- 
tion rt=2.6±0.28 and spiral termination radius rc=5.38±0.54 are 
represented by horizontal error bars at the top. 

models of flp (r) , we express the a priori assumption that 
the bar pattern speed is constant by using an S hke that 
in equation (|14p but where, for n > t, S is reflective of 
either a constant, linear, or quadratically varying spiral 
pattern speed. Solutions with spirals of order 0, 1, and 
2, then, have a total of 4, 5, and 6 degrees of freedom, re- 
spectively. The free parameter t we restrict for all models 
such that 1.8<rt<3.0, according to bar-length estimates 
from Figure[3] (aB~2.2) and Figure ID 

For Ar=0.3 bins, we require a total of 71 slices 
(35 on each side) to reach the edge of the surface 
brightness at r'^10.5. This places the cut bin between 
15<c<20, according to the previously motivated restric- 
tion 4.5<rc<6.0 estimated from Figures [3] and SI 

In light of the discussion in § 13.31 we construct er- 
rors for our estimates to reflect the expected accuracy of 
TWR bar and spiral solutions given a particular observa- 
tional scenario. Specifically, we perform the calculation 
for a range of SAs spanning the upper half-plane of the 
galaxy (quadrants I and II), namely ±15°, ±45°, and 
±75°. Each SA corresponds to a unique disk PA. The 
resultant bar and spiral estimates generated using the 
Ar=0.3 bin width (to be discussed at length in the fol- 
lowing sections) are listed in Table [1] and the average and 
rms of the best-fit solutions for this SA range are shown 
in Figure [5) There, horizontal error bars represent the 
dispersion in rt and Tc in the solutions based on varia- 
tions from SA to SA for the Ar=0.3 radial bin width. 

Rewardingly, the best-fit solutions are quite accurate; 
the comprehensive spiral and bar estimates in Figure [6] 
are 6.7% and 8.3% from their actual values (with ~8% 
error in fli, but a slightly larger error of ~18% in ilg). 
So, too, are the determinations of rt and rc, according to 
Figures [4] and [5] Furthermore, our solutions correctly re- 
produce the functional form of the spiral pattern speed. 
At all SAs, out of all solutions with rt and Vc within their 
restricted ranges, the lowest reduced solution corre- 
sponds to a constant spiral. Figure [7] shows a compari- 
son between the actual and best-reproduced <i;> at each 
slice position yi used in the calculation for the 15° slice 
angle. The <v> reproduced by the best-fit constant bar 
and constant spiral solutions are shown in the left panel, 
while those from the optimum (lowest-x^) solution with 
a quadratic spiral are plotted in the right. 



TABLE 1 

TWR ESTIMATES FOR SIMULATION I. 



SA 


Ui 


ns 


n 




75" 


0.325 


0.149 


3.0 


5.7 


45° 


0.327 


0.111 


2.4 


4.8 


15° 


0.343 


0.184 


2.4 


4.5 


-15° 


0.269 


0.172 


2.7 


4.8 


-45° 


0.313 


0.199 


3.0 


5.4 


-75° 


0.303 


0.203 


2.4 


6.0 




0.29 


0.18 


2.5 


5.0 



Note. — TWR bar and spiral pattern speeds from the barred 
spiral simulation listed here are estimated with TWR solutions 
calculated using a Ar=0.3 bin width for a range of SAs. The third 
and fourth columns list the connate estimates for rt and Tc ■ Values 
for the actual pattern speeds are shown in the last row. 

The close reproduction of the actual pattern speeds by 
the solutions in Figure [6] occasions further evaluation of 
strictly unregularized TWR calculations. Consider the 
values in radial bins inside r'^2.0 in typical unregular- 
ized solutions for this simulation (Figure [2]). That there 
is little to no indication of ilb~0.29 in the left plot is 
perhaps not surprising: the large number of bins which 
accompany the choice of the small bin width would seem 
to guarantee a high level of noise propagated through- 
out the solution. However, in the slightly more stable 
solution with the wider bin width, the inner bins are still 
unrepresentative of the actual pattern speed in this zone. 
We can understand this as a systematic error introduced 
by the noise which not only propagates but also com- 
pounds as the full solution assembles from the outermost 
to the innermost radius; rip(r) in a given bin reflects er- 
rors from all exterior radial bins making the value in that 
bin more likely far removed from the actual value. 

4.2.1. Morphology-dependent Effects and Intrinsic 
Limitations 

The use of simulations which can be studied at multiple 
projections provides us with perhaps the most critical 
assessment for the accuracy of TWR solutions. Figure 
[6] suggests that the TWR method should perform well 
for any given viewing angle. However, though still quite 
small, the rms in each estimate is largely reflective of the 
non-trivial effect that the orientation of the pattern with 
respect to the slice angle used in the calculation can have 
on solutions. 

We can understand the origin of the differences in so- 
lutions for the range of SAs in Table [T] by considering the 
impact of the limited azimuthal range of the bright spiral 
enhancement (clear from the surface brightness distribu- 
tion in Figure [3] where the spiral extends almost per- 
pendicular to the bar major axis). That is, at all slice 
angles the quadrature accumulates fragmentary informa- 
tion from the spiral since only some of the slices that 
cross the full radial zone of the spiral pattern intersect 
the strong spiral structure. But whereas both the bar 
and spiral estimates seem to suffer at SAs in quadrant II 
(i.e. positive SAs), our measurements of Qg in quadrant 
I are quite accurate. According to the morphology, in 
quadrant II it appears that the limited sampling of the 
spiral asymmetry implicit in slices other than those that 
also pass through the bar entails slightly less accurate 
spiral estimation. 

To interpret the distinction between solutions from the 
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Fig. 7. — Comparison of model solution-reproduced (closed circles) to actual (open circles) integrals <v>i=bi/ J T,dx as a function of 
slice position y for the barred spiral simulation for a), the best-fit constant Qf,, constant Qg solution, with rt=2A and rc=4.5 and b). 
constant Qi,, quadratic Qs solution with rt=2.7 and rc=5A for SA=15°. Only those slices which show a contribution from bins inward of 
Tc are shown. The adopted global error (t<"> is shown in the upper right. 



two quadrants, consider the combined influence of regu- 
larization and our chosen weighting scheme. Specifically, 
since the regularizing S induces the coupling of all bins 
within the same radial zone, even when the spiral-zone 
crossing slices do not intersect the strong spiral enhance- 
ment, the bins there inherit information from bins at the 
same radius from slices which do intersect the arms. Ac- 
cording to our weighting scheme, however, this coupling 
is not uniform; the degree of support at each azimuth 
is influenced by the measurement errors for each slice 
which grow larger from out to in. As a result, Qg is best 
constrained by information from slices passing through 
the outer radial zone of the spiral alone (not those pass- 
ing through the bar). So when in quadrant II the spiral 
asymmetry appears in only the inner slices, fig is less 
precise then when these outer slices clearly intersect the 
spiral arms. According to equation ([5]), since the bar es- 
timate is directly related to that of the spiral, the result 
for these SAs is error in both 17^ and fib. 

The corresponding determination for the radial domain 
of the bar pattern, on the other hand, is not as obviously 
sensitive to this issue. Not only is the bar-end reasonably 
well defined in both quadrants, but information from the 
bar which contributes to the parameterization of rt is 
reinforced with regularization in the manner described 
above. The spiral termination radius, too, seems fairly 
consistent from SA to SA. But since the asymmetry is 
weaker (and there is less information) at that location in 
the disk, we find that this parameter requires the most 
restriction (indeed, our determinations of Vc completely 
span the allowed range). 

Overall, then, our determinations for rt and Tc are sta- 
ble and accurate, each with less than 11% error, even 
in quadrant II. Nevertheless, since the pattern speed es- 
timates from SAs in quadrant I seem to comparatively 
benefit from the high quality of information from both 
patterns, we conclude that position angles which provide 
the most uniform slice coverage of all patterns in the disk 
are preferred. 

The largest disparity between step-model solutions 
from various slice angles can largely be attributed to lim- 
itations in determining the location of the transition be- 
tween the two patterns. That is, solutions are affected by 
the finite bin width inherent to the numerical calculation; 



slight incompatibility between the actual transition and 
that to which the solution is limited (given the bin width) 
can result in errors around 10%. This is a more pervasive 
effect than morphology alone and more obvious with the 
use of a slightly larger bin width than Ar=0.3. Nonethe- 
less, we can make several informed inferences about the 
result of the finite bin width by considering the nature of 
the TWR calculation. Specifically, since the transition 
determines the contributions of inner and outer patterns 
to the integral on the right hand side of equations ^ 
and ([7]) through the matrix elements Kij, a mismatch 
between the transition bin and the actual transition ra- 
dius will corrupt the separation of the contributions of 
the two patterns. Note that this is precisely the source of 
the covariance between model parameters intimated in § 

To understand the effects of (minor) radial pattern mis- 
assignment, consider a step-function solution that pa- 
rameterizes the extent of a bar pattern speed along with 
that of a spiral. When the transition bin between the 
two patterns underestimates the actual transition by a 
fraction of the bin width, for instance, we would expect 
the numerical calculation to effectively subtract off from 
bi a contribution from the (lower) outer pattern where a 
higher pattern actually exists (e.g. eq. [8]), slightly rais- 
ing the value of the inner pattern speed. Conversely, we 
expect an overestimation of the transition to result in a 
slightly lower inner pattern speed. More subtle effects 
can occur, however, depending on the geometry of the 
patterns, as illustrated for this simulation in § 14.2.41 

For a range of SAs, we emphasize that the resultant 
errors-in possible combination with an undersamplcd 
transition-are minimal, as long as the bin size is suffi- 
ciently small, and are not the result of vastly different 
transitions for each SA; we find that the transition be- 
tween patterns in the best solutions is generally rela- 
tively stable with changes in slice angle. This may be 
unexpected from the perspective of the traditional TW 
method since the projected length presented to the slices 
by the bar will depend on the relative orientation of the 
slices with its major axis (and slice-orientation errors 
tend to be large, as detailed next for this simulation). 

4.2.2. TWR vs. TW 
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TABLE 3 

Traditional TW estimates for Simulation I. 



-0.2 0.2 

<x> 



Fig. 8. — Plot of <i;> vs. <x> for all slices with \y\ <5.1 for the 
barred spiral simulation at SA=-45°. The dashed line is the best- 
fit straight line to the inner eleven slices (stars), while the dotted 
line is the best- fit straight line to all slices shown (stars and x's). 



TABLE 2 

Optimal TW bar estimates for Simulation I. 



SA 


A' slices 


1 Ip 


a 


75" 


7 


0.270 


±0.010 


45° 


11 


0.282 


±0.062 


15° 


17 


0.208 


±0.007 


-15° 


15 


0.248 


±0.026 


-45° 


11 


0.294 


±0.008 


-75° 


9 


0.233 


±0.053 



Note. — All entries originate through the use of an optimal 
number of slices spaced at Aiy=0.3. The number of slices A'^ used 
in the TW calculation are indicated. 



In light of the above results, we next examine the im- 
provements available to pattern speed estimation using 
the TWR calculation relative to the TW method. We 
specifically compare the bar pattern speed estimates ar- 
rived at using the TWR method with those using fully ex- 
tended TW integrals. Though TW estimates of fif, in the 
presence of a secondary structure may also be attempted 
using truncated integrals that extend to just past the end 
of the bar (such that information from the bar alone is 
dominant), we examine the former case for two reasons: 
1) to compare the two methods under identical conditions 
(i.e. using the same data points <v>i) and 2) to study 
the influence of the relatively weak spiral (and evaluate 
the assumption of negligible non-axisymmetric motions 
beyond the bar). Though we use fully extended inte- 
grals to perform the TW calculation-even when making 
estimates of the bar pattern speed-the innermost slices 
clearly supply evidence for a bar pattern speed that is 
distinct from that of the other structure in the disk. Fig- 
ure [8] shows a typical plot of <v> vs. <a;> for this 
simulation where the inner eleven slices are indeed best 
fit by a steeper slope than for all slices. 

For all of the other SAs studied in the previous section, 
we measure fip with a unique number of bar-crossing 
slices. This is intended to reproduce an optimal TW 
observing strategy that makes use of only those slices 
that intersect the enhanced bar surface density. That is, 
since the projected length presented to the slices by the 
strong bar structure depends on the relative orientation 
of slices with the major axis of the bar, the number of 
slices that intersect the bar enhancement varies from SA 
to SA (from seven to 17 for this simulation at the six 



SA 


inner 


17 slices 


all slices 






cr 


i Ip 


fj 


75" 


0.234 


±0.004 


0.115 


±0.008 


45° 


0.217 


±0.045 


0.094 


±0.018 


15° 


0.208 


±0.007 


0.157 


±0.009 


-15° 


0.216 


±0.016 


0.121 


±0.007 


-45° 


0.278 


±0.009 


0.166 


±0.009 


-75° 


0.175 


±0.010 


0.128 


±0.010 



Note. — As in Tabled the TW estimates for this SA range are 
generated with slices spaced at Aj/=0.3. Here, the second column 
lists the estimates Hp along with errors a using the inner seventeen 
slices, respectively, while the last column is from a fit to all slices 
with |3/1<5.0, out to the inferred spiral-end. 



studied SAs) . Table [2] lists this optimal number of slices 
N at each SA along with the corresponding pattern speed 
estimate. All entries in the table correspond to slopes 
of best-fit straight lines (and the corresponding intrinsic 
scatter) in plots of <v> vs. <x>. Upon inspection. 
Table [3] seems to suggest that even in the presence of the 
spiral asymmetry information from the bar is maximal in 
the bar-crossing TW integrals. Despite also reflecting a 
contribution from the spiral pattern, these bar estimates 
are fairly accurate (though presumably not as accurate 
as would be the case for a strictly SBO galaxy such as 
NGC 7079). However, as in Table [31 if we extend the 
slice coverage at each SA out to |2/|~2.4 using the inner 
17 slices (at spacing) -closer to the full extent of 

the radial zone of the bar pattern, according to Figure[51- 
then the quality of the TW estimates diminishes with the 
inner, bar pattern speed estimates approaching that from 
a fit to all slices with \y\ <5.1; even the true bar-crossing 
slices contain non-negligible information from the spiral 
pattern. 

Unlike the TW method, the TWR method is not rele- 
gated to the use of only those slices where the bar con- 
tribution is maximized. In principle, the inner solutions 
are accessible through the very use of information from 
beyond the bar (namely from the zone of the spiral), and 
are improved when this information is radially coupled 
(e.g. by regularization). This aspect of the calculation 
allows for the return of pattern speed solutions without 
reference to an assumed pattern extent, and provides, 
moreover, an independent means of determining the ra- 
dial domain of patterns. 

Of course, this is not to imply that a deficiency of us- 
able slices in TW estimates prevents accurate pattern 
speed measurement, or determination of pattern extent, 
for that matter; using slices that cross primarily through 
the enhanced emission from the bar, TW estimates from 
SBO galaxies have been successfully used to observation- 
ally confirm that bars end at or inside their corotation 
radii. However, it does suggest that in the presence of 
non-negligible asymmetry exterior to the bar, TW bar 
estimates are susceptible to errors introduced by the use 
of slices positioned near the bar end which are presumed 
to reflect the bar pattern speed but in reality include a 
significant contribution from this structure. Consider a 
typical plot of vs. <a;> for bar-crossing slices. The 
best-fit slope determination for fip from such a plot is 
primarily governed by slices with the largest <f > and 
<x>. Since in SBO galaxies <v> and <a;> approach 
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Fig. 9. — Variation of <x> (top), <v> (center) and 
f2p '^ = <i)>/<a;> (bottom) with for the sHce at y=1.2 with 
SA=15° for the barred spiral simulation. 

zero in slices at or near the projected bar end (since they 
are presumably too far past the strong bar structure to 
contain information about the bar and mark, rather, a 
return to axisymmetry in the disk) , even when TW esti- 
mates consider these 'near zero' slices, they contribute 
minimally to estimates for fip. (See iDebattistal |2003| 
where studies with a simulated SBO include a number 
of slices sampling the full extent of the bar.) But when 
there is considerable asymmetric structure present be- 
yond the end of the bar, similarly positioned slices will 
reflect this contribution, impairing measurement of the 
true bar pattern speed. 

In the barred spiral simulation, this consequence can 
be characterized upon inspection of <x>, <v> and 
^p^—<v>/<x> as a function of the limit of integra- 
tion Xmax^Xo in a typical slice at y=1.2 (see Figure 
[9] for SA=-75°). for this slice, which according 

to our estimate for as crosses the outer region of the 
bar, actually seems more reflective of the spiral pattern 
speed (beyond the discontinuity at Xo~3.2). Indeed, 
past where <x> crosses zero we can infer that 

this bar-crossing slice contains substantial participation 
from the spiral; both <x> and <v> decrease to a dip be- 
tween Xo=4:.0 and 6.0 before reaching a plateau. This is 
evidence that, despite the relative weakness and limited 
extent of the spiral as indicated by the surface brightness 
distribution, the asymmetry in the disk should not be 
considered dominated by the bar alone. Moreover, this is 
a clear indication that the TWR calculation, which iden- 
tifies and effectively removes the spiral contribution from 
the fully extended E-weighted velocity integrals, can im- 
prove upon traditional TW bar estimates that also use 
fully extended integrals. Of course, as displayed by Fig- 
ure ini integrating only between ±Xo~ 2.5 may alone 
provide a reasonable bar pattern speed and relieve all 



other TW slices from the spiral contribution. However, 
this would provide only a single pattern speed estimate 
where two are possible; using the TWR method and fully 
extended integrals, both ilf, and fl^ can be measured. 

Figure [9] also raises a crucial point related to the re- 
quired limit of integration along each slice in the TWR 
calculation: if there is a clear plateau in the integrals 
reached before the edge of the map boundary, why not 
simply truncate the integrals where they have converged 
(common to TW estimates) rather than use integrals 
extending to the edge of the surface brightness and 
which present outermost bins that we demand must be 
'cut' (that is, calculated without regularization) , any- 
way? The plateau reached at Xq^S.O would suggest 
that truncating the integral there could suitably ac- 
count for information from the major sources of asym- 
metry in the disk. However, this same distinction is not 
clearly shared by all shces, especially those at large \y\. 
If we associate the plateau in this slice at j/=1.2 with 
the limit of integration Xmax=&-'2 and hence the total 
radial extent Rmax required of the quadrature where 
Rmax — \J ^max + {v I cosQ;)^, then this locates the out- 
ermost slice position at ymax = {Rmax — I^t) j cos a as well 
as the extent of the integral along this slice. It is easy 
to check that (at least for this slice orientation) the inte- 
gral has not achieved convergence by this point, nor have 
most other integrals in the disk by Rmax- 

Rather than risk ignoring information critical for char- 
acterizing the patterns uniformly throughout the disk, 
then, we choose to include in the quadrature all informa- 
tion out to the edge of the surface density. Indeed, this 
serves to perform a function similar to truncating TW in- 
tegrals. The difference in the two procedures arises from 
the fact that while the limits of integration for each indi- 
vidual slice can be adjusted for a given structure of inter- 
est in the TW calculation, the quadrature in the TWR 
method delineates specific bounds which must encom- 
pass complete information from all extended patterns in 
the disk. 

4.2.3. Systematic Errors 

In this section we use the barred spiral simulation to 
consider the errors introduced to TWR pattern speed 
estimates in real galaxies. We can expect errors in the 
assumed PA^isc to dominate errors in the TWR calcula- 
tion, given that such errors translate significantly to in- 
accuracies in the traditional TW calculation via the line- 
of-sight velocity integral which is also, of course, a promi- 
nent feature of the radial TW equation. Table [3] summa- 
rizes the results for a standard PA error of SpA = ±2° 
on the SAs in quadrant I chosen for their advantages, as 
evidenced by the discussion in § 14.2.11 The average and 
rms for bar and spiral estimates from both the TWR and 
the traditional TW methods are listed. (The TW 'bar' 
estimates are obtained by fitting to the inner, nominal 
number of slices listed in Table [21 while all slices with 
|?/|<5.1 are considered in the 'spiral' estimates.) 

Even this small 6pA introduces considerable errors (rel- 
ative to the known pattern speeds) to both types of bar 
estimates. These errors in fif, can be many times larger 
than the formal rms. But whereas the errors are com- 
parable in the TWR and TW bar estimates, the errors 
in the spiral estimates tend to be smaller with TWR 
than TW. Rewardingly, with this error not only are the 
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TABLE 4 

PA ERRORS IN TW AND TWR ESTIMATES FOR SIMULATION I. 





&PA 


= +2" 


5p,4 


= -2* 




fit 




fit 


fls 


TW* 


0.179 


0.162 


0.354 


0.115 




±0.055 


±0.058 


±0.04 


±0.02 


TWR 


0.214 


0.211 


0.380 


0.158 




±0.031 


±0.045 


±0.042 


±0.023 



Note. — Entries correspond to average bar and spiral estimates 
from SAs in quadrant I (-15°, -45°, and -75°) with PA errors SpA = 
+2° and -2°. Estimates from both the traditional TW (*bar 
estimates using the nominal number of slices for each SA listed 
in Table [2] and TWR methods are listed). 

TWR spiral solutions still definitively constant and ac- 
curate to '^15%, but the radial domain of both pattern 
speeds are still well-determined. The transition between 
the bar and spiral rt and the termination radius of the 
strong spiral pattern rc are effectively unchanged from 
the SpA = case; for both 6pa = +2° and Spa = —2° 
we find rt=2.8±0.37and rc=5.4±0.46. 

Besides the effects on PAiUsc measurement as studied 
by iDebattistal (|2003). galaxy inclination and ellipticity 
play perhaps more prominent roles as sources of error 
in TWR solutions relative to traditional TW estimates. 
Presumably, large inclination errors will prevent the as- 
sociation of information into accurate radial bins, given 
that r = Uobs/ cos a. We expect this effect to be minimal 
at moderate inclinations since dr oc da sin a, and most 
significant at small inclinations where one would gener- 
ally find that the difficulty in inferring in-plane morphol- 
ogy and kinematics makes the the TW method imprac- 
tical in any case. At a moderate, 45° inclination we find 
that the barred spiral solutions for ± 3° inclination error 
differ from the actual pattern speeds by only the change 
in sin a introduced by the line-of-sight velocity. 

4.2.4. Transition Mis-identification 

We have shown that the regularized TWR method can 
be used to parameterize the number and radial domain 
of multiple pattern speeds in a single disk. Formally, the 
contribution of each to the line-of-sight velocity integral 
is established through the designation of a transition be- 
tween patterns. In our scheme this transition is a free 
parameter, but the method, of course, could plausibly 
assimilate other transition-iden tification me thods to sim- 
ilar effect, much like that in the lMacieiews ki (2006') adap- 
tation of the TW method. Specifically, in .Macieiewskil 
()200(][ ) a plateau in the integrals J^Xo ^^^^ with varia- 
tion in Xq is associated with the transition from an inner 
to an outer pattern (see lMacieiewskil ()200Q) for details). 
This transition is then used to separate the disk surface 
brightness into two unique components (one for an inner 
secondary bar, one for an outer primary bar), thereby 
governing the decoupling of the pattern speeds. 

We here pursue this type of diagnostic for the case of 
the simulated barred spiral in order to test the reliabil- 
ity of employing the TWR method with such indepen- 
dent evidence for pattern extent. Figure [TU] plots the 

values of J^Xo ^^^^ / I-Xo '^'^^ ^ function of Rq for 
five bar-crossing slices (7/1=0.0, y2=-0.6, y3=-0.3, y4—0.6, 
y5=0.9) at SA=45° where Ro=^/yf+X§. 
For this projection, there seems to be a plateau at 
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Fig. 10.— Variation of <x>=f^°^T:xdx/ f^^^Edx with 

Ro=^y^ + for five slices at SA=45° (yi=0.0 (dash-dot); y2=- 

0.6 (dash); j/3=-0.3 (soUd); y4=0.6 (dot); j/5=0.9 (dash-dot-dot)) 
in the barred spiral simulation. The plateau range indicates a tran- 
sition between Rq ~1.0-1.5. 



TABLE 5 

TWR ESTIMATES FOR SIMULATION I WITH A MIS-IDENTIFIED 
TRANSITION. 

SA fit fls rt r'c 

076 0093 O 

45° 0.348 0.152 0.9 4.8 

15° 0.367 0.191 1.2 4.5 

-15° 0.336 0.186 1.2 4.5 

-45° 0.336 0.212 1.2 5.7 

-75° 0.313 0.200 1.2 5.7 

Note. — All bar and spiral pattern speeds listed here are esti- 
mated from TWR solutions where rt is restricted to between 0.9 
and 1.5. As in Table[T] solutions are calculated using a Ar=0.3 bin 
width. The third and fourth columns list the connate estimates for 
rt and rc. 



Rq ^ 1.2. We note that this value is smaller than the bar 
pattern extent indicated by our best-fit solutions (and 
the major axis bar length estimated by inspection of 
the surface brightness distribution) and furthermore, the 
same analysis performed over the range of SAs does not 
always as clearly show the same behavior. Presumably, 
this particular value is more indicative of the bar minor 
axis length than the full radial zone of the bar pattern, 
since slices at a 45° SA sample along X perpendicular 
to the bar major axis; other slice orientations are simi- 
larly limited to sampling the bar according to its specific 
projection. 

Allowing that information from the full radial zone of 
the bar is not manifest in this type of indicator, we here 
proceed to assess the consequences for inner and outer 
pattern speed estimates when the transition between the 
two patterns is misidentified. Table [5] lists the TWR bar 
and spiral estimates for each of the SAs studied in the 
previous sections where we have limited the transition to 
0.9<rf<1.5. At all SAs, this error of several bins in rt 
causes an overestimation of the inner pattern speed. The 
outer pattern speed, on the other hand, though slightly 
raised, is still reassuringly accurate. We can interpret 
the inaccuracy in Qi,, then, as the result of the mis- 
association of information from one pattern to the other 
via equation ([8]) , as discussed at the end of § 14.2.11 Addi- 
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Fig. 11. — Face-on display of the spiral simulation's surface 
brightness distribution projected with a -30° rotation about the z 
axis, with orientation as in Figure [3] 




Fig. 12. — Face-on display of the spiral simulation's velocity field 
projected with no rotation about the z axis shown here with the 
kinematical major axis running from left to right. Contours are 
spaced at A«=0.1. 

tionally, we can attribute the greater inaccuracy in solu- 
tions from quadrant II to the reasons discussed in 14. 2. II 
Though the subtleties in Table [5] are most likely specific 
to this simulation, we emphasize that systematic pattern 
speed errors introduced by transition mis-identification 
are generic to the nature of the calculation. 

One of the greatest strengths of the TWR calculation 
is that the transition is in principle a free parameter 
(within limits) and need not be restricted to a single, 
pre-determined value. We therefore recommend letting 
the results of the TWR calculation speak for themselves: 
given sufficient resolution and reasonable measurement 
errors, step-model solutions with the most realistic tran- 
sition should be recognizable by how well they reproduce 
the actual <v>. Since the transition determines the sep- 
aration of the patterns by interpreting the contribution 
made to these integrals by each, the natural result is the 
most accurate determination of the pattern speeds pos- 
sible. 

4.3. Simulation II: Spiral Galaxy 



In the previous section we showed that the TWR 
method is capable of detecting and measuring a constant 
spiral pattern speed that spans less than a third of the 
disk. Here, we test the aptitude of the TWR method 
in measuring a radially varying spiral pattern speed that 
subsists over a large radial zone. Since the strong spiral 
surface density enhancement in this simulation (Figure 
lll[) has only moderate azimuthal range like the spiral 
in § 14. 2[ we further explore the likely limitations intrin- 
sic to detecting spiral nonaxisymmetry with a given slice 
orientation. 

The two-armed spiral featured in this simulation ex- 
tends over a large portion of the disk and is strong 
both in the surface brightness distribution (Figure [TT|) 
and as traced by departures from axisymmetric rotation 
(streaming motions) in the velocity field (Figure [T^ . We 
estimate the extent of the spiral structure from that of 
the dominant m=2 component in the Fourier power spec- 
trum plotted in Figure [T31 With the expectation, then, 
that the spiral structure exists between 0.5<r<3.5, we 
restrict our spiral pattern speed solutions to the radial 
zone rt<r<rc bordered at the innermost and outermost 
radii by two independent sets of unregularized bins. (For 
our models, 0<rt<1.0 and 2.8<rc<4.2.) Between the 
radii marked by the free (though restricted) parameters 
rt and Tc, we allow each spiral pattern speed solution to 
vary with radius as an n*'' order polynomial where n—0 
to 2. 

The two sections of 'place-holding' unregularized bins 
serve to isolate the solution in the radial range of interest; 
by minimizing errors due to either incorrectly incorpo- 
rating or imposing an ill-prescribed form to bins outside 
the radial zone of the spiral, we reduce the introduc- 
tion of inaccuracies to the spiral pattern speed solution. 
The first, outermost section covers the radial zone iden- 
tified by visual inspection of the surface brightness (and 
substantiated by the Fourier power spectrum) where the 
dominant spiral ceases to extend and where the weakness 
of the Fourier components suggests a region without a 
noticeable pattern. As a result, we consider this zone un- 
suitable for reliable extraction of a pattern speed. (This 
last point can also be evidenced by traditional TW values 
from slices that pass solely through the outermost disk; 
the low counts there lead to a large degree of variation in 
Qp estimates from slice to slice). The second, innermost 
section covers the radial zone inside rt~1.0. Within this 
radius we do not expect to be able to extract a realistic 
pattern speed solution since neither the morphology, ve- 
locity field, nor the power spectrum indicate a departure 
from axisymmetry. 

Lacking such evidence in the innermost radii, it would 
be entirely possible to proceed without the use of un- 
regularized inner bins (as would be the case when the 
indicators for such a measure are perhaps less obvious). 
Indeed, allowing an independent regularized solution to 
exist at r<rt, the TWR calculation still performs well; 
since these bins cover a rather small region of the disk 
and thereby contribute minimally to the through only 
the innermost slices, we find that the spiral solution out- 
side of this region closely resembles that in the case where 
the inner bins are unregularized. However, we proceed in 
the manner described above with the expectation that, 
if only minorly, our spiral pattern speed solution will be 
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Fig. 13. — Fourier power spectrum of the spira l simulation's 
surface brightness distribution shown in Figure 1111 Modes up to 
m = 4 are plotted as a function of radius with lines for m = 1 in 
dot, m = 2 in solid, m = 3 in dash-dot-dot, and m = 4 in dash. 
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Fig. 14. — Plot of Hp as a function of radius for the spiral simu- 
lation as derived from the time evolution of the m = 2 component. 
Curves for f2 and f2 ± k/2 arc shown in gray. 

improved. 

The actual pattern speed at the time of the snapshot 
shown in Figure [TT] as derived from the time evolution of 
the phase of the m—l component is plotted as a function 
of radius in Figure 1141 This plot confirms that within 
r~1.0, the pattern speed is ill-quantified, with the values 
for fip at the inner-most radii oscillating between posi- 
tive and negative values outside of the vertical range of 
the plot. At the largest radii, the pattern speed is char- 
acterized by scatter presumably reflective of the lack of 
a noticeable pattern in Figure 1111 

The pattern speed between 1.0<r<3.0 to be repro- 
duced by our solutions shows high-order variation with 
radius. Inside of r~2.0 where the pattern speed is at 
a maximum (rip,maa;~0-3), the pattern seems to be un- 
winding, while at larger radii the pattern speed decreases 
with increasing radius. 

This behavior is only modestly indicated by traditional 
TW estimates. In plots of <v> vs. <x> (see Figure ITSl 
for a comparison of plots generated at four SAs), adjacent 
slices trace out a figure-of-eight shape characteristic of 
complex radial behavior (namely winding). However, the 
evidence for a variable best-fit slope (expected for a pat- 
tern speed that unwinds and winds) is not comprehensive 
or even readily apparent in all cases. (In fact, all slices 
are seemingly well fit with a single slope.) The radial 
dependence of Op becomes more apparent upon inspec- 
tion of the variation of fip with slice position (see Figure 



fT6)) . but again, the radial dependence gets smoothed out 
since fip for each slice is the result of averaging over all 
sampled radii. Furthermore, for the four SAs shown in 
Figure I16i there is no single radial behavior implied by 
all. 

The TWR solutions, on the other hand, are capable 
of clearly displaying high-order variation in r2p(r). The 
average and rms of the best-fit solutions with Ar=0.2 
bins at six diflPcrcnt SAs (±15°, ±45°, and ±75°) is plot- 
ted in Figure [T7] between the average values of rt and 
Tc. For each SA, we find that the best-fit solution is 
quadratically varying, correctly reproducing the winding 
behavior present in the actual pattern speed. However, 
the parameters vt , r,, and the size and location of ^^^max 
differ slightly in each solution. These differences tend to 
reflect the influence of morphology and slice orientation, 
as in § 14.2.11 The value of the peak in rip(r), for ex- 
ample, is higher and more pronounced in solutions from 
quadrant II than quadrant I; at positive (negative) SAs 
^^p,ma£c'^0.34 (0.27), on average. Together with differ- 
ences in the location where d£lp(r)ldr = (which varies 
for the six SAs within the range 1.6<r<2.2), upon aver- 
aging, the result is a slightly under-measured peak value 
occurring at r~1.8 (in accord with the actual location). 

The vertical error bars on bins at large and small r sim- 
ilarly reflect variations in the location of and rt with 
SA. Principally, the solution suffers contamination from 
bins that perhaps yet contain information from the inner 
or outer axisymmetric zones. But we also find that so- 
lutions in quadrant II tend to decrease from ^p.max to a 
value Tc that is further in by about 10% than for solutions 
using slices at the perpendicular orientation. Neverthe- 
less, we find the average rt=0.8±0.2 and rc=3.2±0.3 to 
be in agreement with the bounds of the spiral pattern 
indicated by the disk surface density and its Fourier de- 
composition. 

Despite these PA-dependent effects, the values im- 
phed by the average solution (Op,max, f^p,mjn|r=o.8, and 
f^p, mm |r=3. 2) are accurate to within 5%, 11% and 10%, 
respectively. That we have correctly reproduced the 
high-order variation of r2p(r), regardless of SA, however, 
is perhaps the most remarkable aspect of the TWR so- 
lutions, even though the detectable variation is only at 
the 30% level. 

Naturally, our solution for rip(r) lends itself to pattern 
winding time estimates. With the average values for the 
maximum and minima implied by our solution we can 
estimate the winding time of the pattern according to 

'^wind — 27r/(flp_^(ja; Op, mm)- (-^^) 

For the outer spiral arm, for example, we estimate an 
average time to wind fwind= 71.64 which is less than 10% 
from the actual winding time ru,i„(;=78.53 observed from 
the time evolution of the simulation. (This, of course, 
assumes that Op(r) does not vary over this time.) 

As this simulation would indicate, even without uni- 
form slice coverage, though it may be slightly more dif- 
ficult to determine with confidence the radial domain 
of the pattern (given large errors in rt and Tc), the 
overall shape, or functional form, for f2p(r') can be as- 
certained. Of course, this is largely influenced by the 
adopted measurement errors ct^"^ for each slice and the 
quality of a priori information that can be gathered and 
employed. With larger errors cr^"^ , for instance, the 
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Fig. 15. — Plot of <D> vs. for every other slice spaced at A3/=0.4 for the spiral simulation at four SAs (clockwise from top left: 

45°, 15°, -15°, and -45°). Each slice is labeled by its distance from the galaxy major axis yi. Adjacent apertures in the -45° case are 
connected by a solid line; this SA shows the clearest signature of winding. The dashed line in all plots is the best-fit straight line to all 
apertures shown. 
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Fig. 16. — Plot of r2p=<i'>/<x> as a function of distance from the galaxy major axis y for the spiral simulation at the four SAs in 
Figure [TSl (clockwise from top left: 45°, 15°, -45°, and -15°). Triangles (circles) mark slices on the y>0 (y<0) side of the galaxy. 
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Fig. 17. — The best-fit regularized solution and error bars from 
the TWR method using Ar=0.2 bins as applied to the spiral simu- 
lation averaged over six SAs. The best solution is shown as a series 
of solid lines for each bin with dashed errors. Errors for the tran- 
sitions rt=0.8it0.2 and rc=3.2±0.3 are represented by horizontal 
error bars at the top. 




Fig. 18. — Face-on display of the double barred simulation's 
surface brightness distribution projected without rotation about 
the z axis, highlighting the inner two bars. 

criterion becomes less discriminating, and it may be 
difficult to distinguish between several different radial 
dependences for rip(r). Additionally, without clear ev- 
idence that limits where the spiral pattern terminates, 
we risk niisidentifying intrinsic radial variation. Indeed, 
if we restrict Tc to less than that implied by the best- 
fit solutions-and search instead for solutions at a second 
Tc-x^ minimum-the pattern speed solutions for all six 
SAs are constant between ft'^0.47 and fc~2.47. As may 
be expected, the average value for this constant pattern 
speed 17p=0.236±0.051 is similar to that suggested by 
traditional TW estimates where Op^TH'=0.207±0.046 on 
average. 

4.4. Simulation III: Double Barred Galaxy 

In this section we address the use of the TWR method 
for the purposes of nuclear bar detection and measure- 
ment using the double barred SBO simulation pictured 
in Figure [TSl In performing the regularized TWR cal- 
culation we again act under the assumption of multi- 
ple patterns in distinct radial zones. Our models for 




Fig. 19. — Fourier power spectrum of the double barred simu- 
lation's surface brightness distribution shown in Figure [181 Modes 
up to m = 4 are plotted as a function of radius with lines for m = 1 
in dot, m = 2 in solid, m = 3 in dash-dot-dot, and m = 4 in dash. 

flp{r) parameterize unique, constant pattern speeds for 
both the primary and secondary bars, known to have 
pattern speeds of flpb=0.23 and 0^6=0.41, respectively. 
From inspection of the surface density and its Fourier 
decomposition (for which the power spectrum is plotted 
in Figure fT9|. we associate the drop in power of the m=2 
component at r^O.S and again at r^^S.O with the end of 
each bar. In step-models for f2p(r), then, we restrict the 
secondary-to-primary bar and primary bar-to-disk tran- 
sitions to within 0.4<r(.i<1.2 and 2.5<rt^2<3.2, respec- 
tively. 

To isolate the bars from the rest of the disk we extend 
the quadrature to the edge of the surface density and 
employ our 'cut' procedure in light of the argument set 
forth at the end of i! l4.2.2l This is particularly compelling 
here since, despite the apparent axisymmetry beyond the 
primary bar in this SBO simulation, the Fourier decom- 
position shows power in the m—2 mode beyond r^3.0, 
especially in the last third of the disk. If the asymmetry 
in this radial zone (which appears only very weakly in the 
surface density) is sustained by an ill-defined, non-unique 
or perhaps unrigid pattern-for which the TW/TWR as- 
sumptions break-down-then associating it with a mea- 
surable pattern speed will likely introduce regularization- 
induced bias to the solution for the interior patterns of 
interest. 

Critically, asymmetry such as this may prevent clear 
integral convergence beyond the primary bar end. So, 
too, can its presence in the integrals be expected to de- 
grade the reliability of pattern speed estimates for the 
structures of interest. Removing the influence of the in- 
formation in this outer radial zone by calculating the 
bins there without regularization is our best chance for 
accurate pattern speed measurement. Given that the de- 
parture from axisymmetry manifest by a small nuclear 
bar will be relatively minor compared with that of other 
patterns in the disk, this is especially relevant for accu- 
rate measurement of fls- In this case, non-axisymmetry 
on a comparable scale may easily upset this structure's 
contribution to the integrals and, if prescribed an incor- 
rect pattern speed model, may introduce consequences 
for the inner most bins in solutions. 

Our procedure for this simulation, however, does not 
quite involve calculating without regularization all bins 
up to the patterns of interest (i.e. to the end of the pri- 
mary bar). As in all cases, but particularly here where 
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we are compelled to 'cut' approximately two thirds of the 
disk, leaving a large portion of the total bins unregular- 
ized may begin to reintroduce unamendable propagating 
noise. To reduce this risk, careful attention has been 
paid to the development and testing of models which 
prevent the destabilization of solutions. (Indeed, the ap- 
propriate balance between noise and stability in models 
for ilp(r) must be explored on a galaxy by galaxy basis.) 
Here, we find that the most stable and realistic models 
for rip(r) are those which include a third, constant pat- 
tern past the end of the primary bar. From this wc might 
infer that, though weak and difficult to discern in Fig- 
ure [181 there exists a spiral pattern outside the primary 
bar, perhaps corresponding to the m=2 component be- 
yond r'^3.0 which remains clear, though modest out to 
r'^b.O. Indeed, we assume that the third minimum in 
the power of the m=2 component at this radius corre- 
sponds to the end of the spiral pattern, and moreover, 
since counts are low in the rest of the disk, that the bins 
beyond r~5.0 are best calculated without regularization. 
We note, however, that we do not necessarily expect to 
measure a realistic pattern speed in this third radial zone. 

Compiling this evidence for two bars and a possible 
spiral we search for the best-fit solutions parameterizing 
two constant pattern speeds (one for the primary bar, one 
for the inner secondary bar) out to the end of the primary 
bar, in addition to a third constant pattern speed fig 
restricted to extend out to 4.3<rc<5.3. We find that 
solutions generated in this manner provide much more 
accurate estimates for the primary and secondary bars 
compared with solutions that are either regularized over 
the full extent of the surface density or unregularized up 
to the end of the primary bar. 

For the purposes of further establishing favorable con- 
ditions for nuclear bar detection, we adopt a small bin 
width Ar=0.15 in the quadrature. As described in ii l4.2.1l 
and § 14.2.41 the mismatch by a fraction of a bin width or 
more between the actual transition and that to which the 
solution is confined can have consequences for both inner 
and outer TWR pattern speed estimates. Real nuclear 
bars will need to be well resolved in order to accurately 
separate the contributions of the two bars. 

The secondary bar i n this simulation is know n to be 
non-rigidly rotating : in [Shen fc Debattistal (|2007f) and in 
iDebattista fc ShenI (|2007D . the amplitudes and pattern 
speeds of secondary bars formed in purely coUisionlcss A^- 
body simulations through the introduction of a rotating 
pseudobulge oscillate as the bars rota te through the com- 
panion s in which they are nested. In IShen fc Debattistal 
(|2007f ). TW estimates of the secondary bar pattern speed 
measured using bulge-only kinematics (the bulge sup- 
ports the nuclear bar alone and a primary bar contri- 
bution need not be accounted for in the TW integrals) 
are subject to marked errors consistent with an origin 
in non-rigid rotation. These errors result in estimates of 
flsb too high on one side of the galaxy and too low on 
the other, in accord with being a manifestation of the os- 
cillations driving radial pulsations which contribute with 
different signs on the two sides of the galaxy. Cancella- 
tion between measurements from both sides of the galaxy 
in global regularized solutions reduces the effect of the 
oscillations. In the discussion to follow, our focus is on 
sources of error in the TWR calculation other than this 
intrinsic effect. 



TABLE 6 

TWR ESTIMATES FOR SIMULATION III. 

SA Qsb ^pb n.i rt,2 Tc 

5:420 (1:231 OUB 0?75 3J5 OS" 

45° 0.548 0.278 0.178 0.75 3.15 4.35 

15° 0.407 0.231 0.063 0.6 2.7 4.8 

-15° 0.404 0.264 0.087 0.75 2.7 5.25 

-45° 0.477 0.259 0.233 0.9 3.15 4.35 

-75° 0.375 0.225 0.028 0.9 3.15 4.5 

- 0.41 0.23 - 0.8 3.0 5.0 

Note. — The secondary bar, primary bar and spiral pattern 
speeds listed here are estimated from TWR solutions calculated 
using a Ar=0.15 bin width for a range of SAs. The last three 
columns list the connate estimates for rt^i, rt,2 and Tc- Values for 
the actual pattern speeds are shown in the last row. 
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Fig. 20. — The best-fit solution and error bars from the TWR 
method as applied to the double barred simulation averaged over 
six SAs using Ar=0.15 bins. The secondary bar f2s(,=0.439±0.058, 
primary bar r2pj,=0.248±0.02, and spiral Qs=0.123±0.07 are 
shown as solid lines with dashed errors. Horizontal error bars 
at the top indicate the dispersion in the secondary-to-primary 
bar transition rf i=0.75±0.09, the primary bar-to-spiral transition 
rt 2=3.0it0.21, and spiral termination radius rc=4.6±0.33. 

The TWR estimates from solutions with Ar=0.15 bins 
at six SAs (±15°, ±45°, and ±75°) are listed in Table 
[6] and the average and rms for this SA range is shown 
in Figure 1201 Compared with solutions generated using 
Ar=0.3 bins in the calculation, these solutions are much 
more stable and accurate; at less than 8% from their 
actual values, the comprehensive primary and secondary 
bar pattern speeds evidently benefit from the integrity 
of the calculated transitions between patterns (we find 
rt,i=0.75±0.09, rt,2=3.0±0.21, and rc=4.6±0.33). 

Though the TWR estimates in Figure!^ are accurate, 
the rms in each suggests 8 to 13% error on the part of the 
method. In addition to the intrinsic non-rigid rotation- 
error expected for the secondary bar measurement, this 
may yet indicate the slight misdesignation of transitions 
between patterns, and we can also assume that (though 
we do not necessarily consider its measurement realistic) 
the third measured pattern speed has a non-trivial influ- 
ence on the inner pattern speeds, as must Qpi, on flsb, for 
instance. 

With so many free parameters in our model-more than 
in any other used in this paper-identifying unequivocally 
the dependence of each on the others is difficult (indeed, 
the rms in each estimate likely refiects the combined con- 
sequence of many such inter-dependencies). As in pre- 
vious sections, however, we can begin to understand the 
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largest source of variability in the estimates from SA to 
SA by considering the influence of morphology and slice 
orientation on the solutions. If we consider that the ma- 
jor axes of the two bars are oriented at ~45° from one 
another at this time step (that is, with the secondary 
bar aligned along the x-axis, the primary bar is a ligned 
along the line y=x, or in the IShen fc Debattistal ()2007f ) 
convention, A(/)=45°), then slices at ±45° are oriented 
along the primary bar minor/major axis. In this case, 
even in the presence of the secondary bar <v> and <a;> 
for the inner primary-bar crossing slices fluctuate noisily 
about zero. 

This behavior is more than enough to prevent accurate 
TW es timates for the primary bar fwhich lRand fc WallinI 
((2OOI find are less reliable when bars are aligned within 
±20° of a principle axis), and so it is perhaps not sur- 
prising that our TWR estimates at ±45° SAs show 
rather large inaccuracies. However, what we find re- 
quires slightly more interpretation: it is the secondary 
bar estimate, and not Slph, that is compromised. 

Unlike the results in § 14.21 with which we could ar- 
gue that the regularized TWR method is, in principle, 
fairly equipped to accommodate any alignment of pat- 
terns (and where even slices at ±15° from the bar axis 
return accurate TWR estimates), here we are not only 
dealing with more model parameters (and their intrinsic 
covariance), but the innermost pattern now has a much 
smaller extent than the structures in the disk. Conse- 
quently, the secondary bar occupies the smallest fraction 
of bins in the innermost slices where, moreover, the errors 
Gi are the largest. Presumably, in this exploitable state 
the secondary bar estimate is sacrificed for the primary 
bar fipf, in the regularized calculation; at SA=±45°, our 
measurement of is overestimated by ~25%, while 
ripb is good to within 11%, on average. (We can as- 
sume that such imprecision also reflects the effect of the 
secondary bar's non-rigid rotation.) At other slice ori- 
entations where the primary bar is more favorably sam- 
pled, the trade-off is less severe (though modestly PA- 
dependent) . 

An impromptu calcul ation using particle s in th e bulge 
only (motivated by the IShen fc Debattistal ()2007[ ) strat- 
egy) seems to confirm the influence of the primary bar in 
the TWR calculation at SA=45°. The error in our lone 
estimate r2sb=0.42 in this case is significantly less than 
when the primary bar is present and also comparable to 
that in bulge-only estimates at other SAs; for the three 
SAs in quadrant II tested, we find risb=0.43±0.01. 

Even if the TWR method cannot accurately measure 
the pattern speed of a pulsating nuclear bar-aside from 
with global regularized solutions, in particular-it is nev- 
ertheless appropriate for use in characterizing and con- 
straining ^sb, especially in the presence of a primary 
bar. The comparative worth of TWR solutions can be 
deduced in light of the distinction between the tradi- 
tional TW and the TWR methods. Our TWR pat- 
tern speed measurement for the primary bar is compara- 
ble to that from the traditional TW calculation; using 
slices covering \y\ <2.25 (such that the bar contribu- 
tion is maximal in all slices), the TW method returns 
l]p6 -r^=0.214±0.011. The TWR secondary bar esti- 
mate, on the other hand, greatly improves upon the TW 
estimate f2s6,TW=0.338±0.149 available using slices cov- 
ering \y\ <0.6. 



Since the nuclear bar appears in only a very small num- 
ber of slices and its contribution is, moreover, easily over- 
whelmed in the extended TW integrals, the TW method 
only modestly recovers Vlgb- In the TWR calculation, 
however, the few bins that cover the nuclear bar are sup- 
plied with large weights Kij compared with the rest of 
the information along each slice. Together with a suit- 
able bin width, this effectively isolates information from 
throughout the nuclear bar extent from that of all other 
patterns in the disk. The quality with which ilgb can be 
constrained then depends on how well the contribution 
from these other patterns is identified and removed from 
slices intersecting the nuclear bar, and moreover relies on 
regularization over both sides of the galaxy to construct 
a global solution. 

Aside from the unavoidable error introduced by the 
non-rigid rotation of nuclear bars expected from N-body 
simulations, from our current study it seems likely that 
similar calculations may be limited according to the de- 
gree of resolution and the relative size of the secondary 
bar. That is, the former dictates the precision with which 
the secondary-to-primary bar transition can be deter- 
mined, and the latter sets the leverage supplied by the 
secondary bar to the x^- Furthermore, given that the 
measurement of an inner, nuclear bar pattern speed is 
subsequently largely affected by how well the other pat- 
terns present in the disk can be measured, success with 
the TWR method requires appropriate models for rip(r) 
that account for such concerns. 

5. SUMMARY 

5.1. Caveats in Applications to Real Galaxies 

The successes of the TWR method as applied to sim- 
ulated data should be obtainable with real galaxies, pro- 
vided that adequate attention is paid to several consid- 
erations. Since the calculation foremost requires integra- 
tion along slices parallel to the galaxy major axis which 
reflect information from all patterns of interest in the 
disk, observations must be able to present a number of 
these. As governed by the resolution or sampling of the 
data, the placement of these slices must be able to define 
a quadrature wherein the position y of the outermost slice 
corresponds to that of the last matrix element along the 
slice requiring the largest limit of integration X^ax for 
integral convergence; when convergence can be reached 
only clearly at the map boundary (perhaps far beyond 
the extent of the patterns), this is essential for correctly 
accounting for all information along each slice. 

Establishing the slice positions and orientation for an 
accurate quadrature also clearly requires superior knowl- 
edge of kinematic and morphological parameters. In gen- 
eral, the same restrictions to the quality of input data in 
traditional TW calculations apply to the TWR method. 
As described in § I4.2.3[ errors in the assumed PA which 
impair TW estimates can also introduce considerable in- 
accuracy to the TWR solutions. The inclination angle 
of the disk must also be well- determined and should be 
preferably restricted to only moderate values (which, be- 
yond the observational requirements of the TW method, 
will keep errors from corrupting the association of infor- 
mation into accurate radial bins). 

Additionally, though more easily overcome in tradi- 
tional TW estimates (.Merrifield fc Kuiiken|[l995l ). errors 
in the systemic velocity and galaxy center might prove 
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critical for the TWR solutions since each side of the disk 
provides an independent solution for Slp(r). Care must 
be taken not to impair the prevailing symmetry along 
each slice, given that each integral bears more than one 
estimate. The incorrect placement of radial bins accord- 
ing to a mis-assigned galaxy center position, for exam- 
ple, could significantly over- or under-estimate the actual 
pattern speed, and moreover make assessment of the true 
radial variation unlikely. 

The regularization procedure developed here itself 
makes further demands on the quality and amount of in- 
formation necessary to perform the calculation. But by 
keeping the amount of information beyond that required 
of the TW method to a minimum, and using standard di- 
agnostics such as Fourier decomposition, the requisite set 
of a priori assumptions can be invoked quite reliably. As 
long the information is accessible, requiring at the least 
theoretical motivation to develop testable models, and 
limited in principle only by the quality of information 
from which it is to be gathered, then the regularization 
should proceed without impediment. 

Of course, unlike the simulated galaxies studied here 
where there is plenty of signal throughout the disk, ob- 
servations of real galaxies may present sensitivity issues. 
While regions of low signal-to-noise in the outer disk can 
be superseded using the 'cut' procedure developed here, 
high quality information from the rest of the disk is an 
obvious priority for the method; the departures from ax- 
isymmetry induced by all patterns present in the disk 
must be clearly detectable. Not only does the calculation 
depend on the presence (or lack) of these signatures-in 
both the surface density and in the velocity field-but the 
mere identification of the number and domain of patterns 
is critical for developing appropriate models for il,p{r). 

This latter necessity may be hard met since, for in- 
stance, it will be rare to observe galaxies with surface 
densities that can be Fourier-decomposed as cleanly as 
is possible with simulated galaxies. Furthermore, unlike 
simulations, it is impossible to establish whether or not 
there exists more than a single pattern speed at each 
measured radius in real galaxies. Since the models de- 
veloped with regularization here are incompatible with 
non-unique pattern speeds, for real galaxies, the choice 
of when and where to consider a transition or to keep so- 
lutions unregularized may be based on a more qualitative 
assessment of where clear structure ends. 

Inevitably, the combination of the above considerations 
(related, overall, to the quality of the data) will deter- 
mine the extent to which the model for the true radial be- 
havior of f2p(r) can be differentiated from other models. 
That is, the criterion with which we judge the good- 
ness of solutions becomes less discriminating the larger 
the measurement error cr^"-^. Since the adopted mea- 
surement error (t<^> for each slice used in the calcula- 
tion and in the estimator must necessarily incorporate 
observational errors based on random noise in the data, 
with severe enough errors different model solutions from 
real, imperfect data may be indistinguishable. 

In addition, systematic errors (likely dominated by PA 
uncertainty) will undeniably challenge the accuracy of 
solutions. In all applications of the method it is criti- 
cal to assess the influence of these errors through direct 
tests of the sensitivity of solutions to departures from the 



nominal values of PA, inclination, and kinematic center, 
for instance. Clearly, this makes 2-D coverage desirable; 
here, 2-3° uncertainties in the PA alone are shown to in- 
troduce around 15% error in measurements of ^p(r) for 
the barred spiral simulation. 

Insufficient resolution or sampling may also impair 
TWR solutions from real galaxies. A large adopted bin 
width not only limits the detectable radial variation in 
r2p(r), but also restricts how well multiple patterns can 
be separated in the resultant quadrature; a mismatch 
between the actual transition and that to which the so- 
lution is confined can have consequences for the estimates 
of both inner and outer patterns. Naturally, depending 
on the models to be tested and relative size of the disk, 
a resolution-constrained bin width is not guaranteed to 
impair solutions for all galaxies. We nonetheless fore- 
see that the only true way to preserve the integrity of 
solutions is with high-resolution observations. 

Data cubes lend themselves well to analysis with the 
TWR method, since unlike long-slit spectroscopic obser- 
vations, the galaxy PA, inclination, systemic velocity and 
(kinematic) center can be derived with errors from the 
data using a tilted-ring analysis on the first moment of 
the cube. Additionally, multiple slices can be defined 
with a single observation. So between radio and sub-mm 
investigations of spiral structure, for instance, and IFU 
spectroscopy with which double bar systems (and even- 
tually double barred spirals, given larger, more sensitive 
IFUs) can be studied, applications of the TWR method 
could be extensive. 

Of course, like all applications of the TW method, the 
observed tracer must be found to obey continuity and 
the relation between the intensity of the tracer and the 
surface density must be linear or well-determined every- 
where. Reviews of several possible tracers argued to suit- 
ably obey the TW continuity requirement can be found 
th roughout the literature , but w e note here that the work 
of iGerssen fc Debattistal ()2007[ ) studying the effect of 
dust on TW measurements of bars may find meaningful 
extension in future TWR studies of multiple patterns in 
late-type galaxies. There, model dust lane features asso- 
ciated with bars introduce e rrors on the order of 20-40% 
(jGerssen fc Debattistall2007t ) . In addition to these errors, 
TWR solutions could possibly be prone to increased er- 
ror from spiral dust lanes at larger radii. Though it is 
beyond the scope of this work to make a detailed as- 
sessment of the sensitivity of TWR solutions to dust, we 
argue that such noise could well be mitigated through 
the use of regularization, and expect no greater an effect 
in TWR measurements than TW, which moreover, will 
be apparent with the use of only optical tracers. 

In the immediate future, we plan to apply the TWR 
method to several high-resolution BIMA SONG CO ob- 
servations of molecule dominated galaxies to search for 
spiral winding, relations between bar and spiral pattern 
speeds, and spiral-spiral mode coupling. (These observa- 
tions include single-dish data and therefore do not suffer 
from missing flux which would be a violation of the con- 
tinuity requirement.) For those galaxies with ISMs not 
dominated by molecular gas, we plan to combine the CO 
with HI data to make total column density maps (as- 
suming the ionized component is negligible). Since the 
CO-H2 conversion factor is critical in combining the CO 
and HI maps (and, of course, in establishing molecule 
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dominance in the former case), it will be necessary to 
test the sensitivity of the TWR method to the adopted 
conversion factor for such combinations. 

Additionally, since warped disks (common in HI) are 
a violation of the TW assumptions, we will also perform 
tests to determine if our cut-off scheme can be used to 
circumvent the warp and thereby extract solutions from 
the rest of the disk. For these applications of the TWR 
method, we plan to construct measurement errors a for 
each slice that reflect uncertainties related to the flux 
cutoff chosen in creating the moment maps. The effect 
of PA and other systematic errors will be assessed by 
testing the sensitivity of solutions to departures from the 
nominal values. 

5.2. Applying the TWR Method 

Although in the interest of testing our strategy for each 
simulated galaxy is somewhat tailored to its unique prop- 
erties, with the above caveats in mind our studies have 
enabled us to develop a general and reasonably objective 
prescription for applying the regularized TWR calcula- 
tion: 

1. Establish the bin width and the corresponding 
number of slices (not necessarily uniformly spaced) 
that are required to achieve converged integrals us- 
ing an X iV quadrature. For the purposes of 
measuring multiple patterns in a single disk, this 
will likely extend to the map boundary. 

2. Compile a priori information by inspecting the sur- 
face density, its Fourier decomposition, and the ve- 
locity field for indications of patterns and to es- 
tablish the expected number and domain of mea- 
surable pattern speeds. This should include the 
identification of regions in the disk susceptible to 
regularization- induced bias. 

3. Develop theoretically and/or observationally moti- 
vated models which parameterize flp{r) according 
to the a priori information. 

4. Incorporate measurement errors into a single cr^"^ 
(and a) for each slice. These should represent 
uncertainties in the adopted intensity noise level, 
and/or other random noise-related errors; system- 
atic errors are preferably determined through di- 
rect testing (see item 9). 

5. Develop the weighting scheme for a reduced es- 
timator which accounts for the total degrees of free- 
dom for the models to be tested. This should re- 
flect expectations for which slices, if any, are most 
critically to be reproduced by the models. 

6. Generate preliminary solutions for the models. At 
this point, the degree of regularization required to 
return solutions according to type should be ex- 
plored. 

7. With finahzed solutions, use equation ([7|) to gen- 
erate a complete set of <v> for each and calculate 
the corresponding xt- 

8. Use the xt to identify the best-fit solution. 



9. Test the sensitivity of the results to other sys- 
tematic effects peculiar to the observation, e.g. 
adopted PA and/or the CO-H2 conversion factor, 
for instance. 

6. CONCLUSION 

In this paper we have shown that regularizing the 
TWR calculation is an effective means of smoothing in- 
trinsically noisy solutions for more precise measurement 
of rip(r). Specifically, (barring a large, limiting resolu- 
tion) regularization admits the use of a much smaller bin 
width than that required to achieve comparable smooth- 
ness in the unregularized calculation. This affords im- 
proved assessment of radial variation as well as more ac- 
curate determination of the transitions between multiple 
pattern speeds (and thus of the values of the pattern 
speeds themselves, in principle). Moreover, with the reg- 
ularized TWR calculation, different theoretically and ob- 
servationally motivated models for the radial dependence 
of Qp (r) can be tested in fairly short time and with only 
the minor addition of information compared to the un- 
regularized TWR and TW methods. 

With a simple scheme for generating n*'* order polyno- 
mial solutions which can be incorporated into step mod- 
els, we have shown that the TWR method can be used 
to parameterize the radial domains of multiple pattern 
speeds. Together with a priori information identifying 
zones in the disk which may be incompatible with mea- 
surement (either because they are characterized by low 
signal-to-noise or show no evidence for a pattern) , we can 
further constrain the extent of patterns while optimally 
reducing regularization-induced bias in pattern speed so- 
lutions. 

As applied to three simulated galaxies, we find that 
the TWR method developed in this manner performs 
with a high degree of accuracy (with less than 15% er- 
ror) both in measurement and in extracting information 
about the true functional form for the pattern speed. 
Tests on a simulation of a barred spiral galaxy indicate 
that not only can the constant pattern speed for a rel- 
atively weak spiral be reliably reproduced, but informa- 
tion about both the pattern speed and the radial extent 
of the bar pattern can also be extracted. (Indeed, we find 
that the bar pattern speed estimate is strengthened by 
the proper use of information from beyond the bar end.) 
And though the bar pattern speed estimates is highly 
susceptible to systematic errors-with PA errors introduc- 
ing the largest uncertainty to TWR pattern speeds, as 
with TW estimates-we find that the identification of the 
transition between the two is relatively stable. 

The TWR method can also be effectively employed 
to measure patterns that are winding in nature. In a 
simulation of a two-armed spiral, the best-fit TWR solu- 
tions from several slice orientations correctly reproduce 
the high-order radial variation of the pattern speed, de- 
spite modest indication that not all orientations supply 
the same authority (this, of course, would seem to depend 
on the morphology of this spiral, in particular). Indeed, 
though the TWR method can, in principle, handle any 
(presumably random) alignment of patterns, in all of the 
simulations studied, slice orientations which provide the 
most uniform coverage of the patterns are preferred. This 
is of particular importance for nuclear bar pattern speed 
measurement, as found in tests of the method on a dou- 
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ble barred simulation. Since the innermost bins which 
provide the foremost leverage on the nuclear bar are also 
the most susceptible to errors from throughout the disk, 
confident measurement requires all other patterns to be 
well constrained. 

In principle, comparable accuracy should be achievable 
on real galaxies. However, these tests do not constrain 
how well the TWR method can perform under severe ob- 
servational limitations which may commonly arise. Not 
only can determinations of the PA, inclination, and dy- 
namical center be subject to considerable errors given 
low-q\iality data, but identifying constraints on the pat- 
terns present in the disk to be incorporated into models 
for rip(r) could prove challenging. Additionally, though 
regularization can reduce the impact of noise on solu- 
tions, large measurement errors for each slice could make 
discriminating between several possible models for Op (r) 
difficult. And most critically, since the nature of the nu- 
merical calculation relies on a relatively small bin width 
to achieve its greatest accuracy, without high resolution, 
some observations may not afford practical solutions. 

Nevertheless, if restricted to high resolution obser- 
vations with adequate sensitivity, and given radially 



stable kinematic parameters, TWR solutions can be 
used to study the connection between multiple patterns 
and the nature of spiral winding. So, too, can we 
expect progressively more satisfactory applications of 
the method; though the number of galaxies to which 
the method can be successfully applied is limited by 
the current generation of instruments, in the future, 
larger IFUs, ALMA, and eventually, SKA should yield 
much higher quality data with larger areal coverage and 
higher angular resolution. This prospect in itself should 
warrant future studies with the TWR method. 
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